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Abstract 

We present the first explicit formulas for the complete set of one-loop helicity amplitudes 
necessary for computing next-to-leading order corrections for e + e~ annihilation into four jets, for 
W, Z or Drell-Yan production in association with two jets at hadron colliders, and for three-jet 
production in deeply inelastic scattering experiments. We include a simpler form of the previously 
published amplitudes for e + e~ to four quarks. We obtain the amplitudes using their analytic 
properties to constrain their form. Systematically eliminating spurious poles from the amplitudes 
leads to relatively compact results. 
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1. Introduction 



The discovery of new physics at colliders relies to a large extent on our ability to understand the 
known physics producing the bulk of the data. For processes involving hadronic jets, perturbative 
QCD predictions are required. Leading-order calculations often reproduce the shapes of distribu- 
tions well but suffer from practical and conceptual problems whose resolution requires the use of 
next-to-leading order (NLO) calculations. In many processes at modern colliders, the dominant 
theoretical uncertainties are due to unknown higher-order perturbative corrections. These correc- 
tions can be enhanced by various logarithms. For some processes, NLO corrections are known, 
and programs implementing them have already played an important role in analyzing data from a 
variety of high-energy collider experiments. Other processes have awaited the computation of the 
required one- loop matrix elements. 

The first type of logarithm contributing to theoretical uncertainties is 'ultraviolet' in nature. 
Such logarithms are connected with the renormalization scale (i, which we are forced to intro- 
duce in order to define the running coupling, a s (fi). Physical quantities, such as cross-sections 
or differential cross-sections, should be independent of ft. When we compute such a quantity in 
perturbation theory, however, we necessarily truncate its expansion in a s , and this introduces a 
spurious dependence on \i. Together, these effects can lead to anywhere from a 30% to a factor of 
2-3 normalization uncertainty in predictions of experimentally-measured distributions. In general, 
the inclusion of NLO corrections significantly reduces the over-all sensitivity of a prediction to 
variations in fi. 

The other type of logarithm is 'infrared' in nature. Such logarithms are connected with the 
presence of soft and collinear radiation. Jets in a detector consist of a spray of hadrons spread over 
a finite region of phase space. Experimental measurements of jet distributions depend on resolution 
parameters, such as the jet cone size and minimum transverse energy. In a leading-order calculation, 
jets are modeled by lone partons. As a result, these predictions either lack a dependence on these 
parameters or have an incorrect dependence on them. In addition, the internal structure of a jet 
cannot be predicted at all. 

In the case of e + e~ annihilation into jets, leading-order predictions for the production of up to 
five jets have been available for quite some time [1,2,3,4,5]. The NLO matrix elements for three-jet 
production and other 0(a s ) observables are also known [3], and numerical programs implementing 
these corrections [6,7] have been widely used to extract a precise value of a s from hadronic event 
shapes at the Z pole [8]. 

Next-to-leading order corrections for more complicated processes are important, however, if we 
wish to use QCD to probe for new physics in other standard model processes. In e + e~ annihilation, 
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for example, four-jet production is the lowest-order process in which the quark and gluon color 
charges can be measured independently. Four-jet production is thus sensitive to the presence of 
light colored fermions such as gluinos [9]. At LEP 2 the process e + e~ — ► (j*,Z) — ► 4 jets is a 
background to threshold production of W pairs, when both Ws decay hadronically. 

The calculation of e + e~ — ► 4 jets at NLO requires the tree-level amplitudes for e + e~ — ► 5 par- 
tons [4,5] (at NLO two of the partons may appear inside a single jet), and the one- loop amplitudes 
for e + e~ — > 4 partons. In a previous paper [10] we presented the one-loop helicity amplitudes for 
electron-positron annihilation into four massless quarks, e + e~ — > (7*,Z) — ► qqQQ (q,Q may have 
the same or different flavors). 

In this paper, we present the e + e~ — > (7*, Z) — > 9955 one-loop helicity amplitudes, as well as 
simplified versions of the e + e~ — > (7* , Z) — > qqQQ amplitudes. We give all contributions at order 
<7 4 in the strong coupling, including those where the vector boson couples directly to a quark loop 
via a vector or axial- vector coupling. We take all quarks to be massless, except for the top quark, 
whose virtual effects we include through order 1 / . Thus the list of helicity amplitudes required 
to construct a numerical program for e + e~ — ► 4 jets is now complete. Indeed, the amplitudes 
presented here and in refs. [11,10] have already been incorporated into the first NLO program for 
e + e~ —> 4 jets [12]. Crossing symmetry and simple coupling constant modifications allow the 
same amplitudes to be used in NLO computations of the production of a vector boson (W, Z, 
or Drell-Yan pair) in association with two jets at hadron colliders, and three-jet production at ep 
colliders. Finally, these amplitudes will also enter the next-to-next-to-leading (NNLO) study of 
three-jet production at the Z pole. Such a study (which awaits the computation of appropriate 
two-loop matrix elements as well) would be desirable in order to reduce the theoretical uncertainties 
in determining a s via this process. 

Glover and Miller [13] have reported on a calculation of the helicity-summed interference 
term between four-quark one-loop matrix elements and the appropriate tree-level matrix elements. 
Recently, Campbell, Glover and Miller [14] have also calculated the analogous 'squared' matrix 
elements for the two-quark two-gluon final state. In neither of these papers did the authors provide 
any explicit formulae. In order to compare with the results described in refs. [13,14], we considered 
the case of virtual photon exchange, dropped the contributions where the photon couples directly 
to a quark loop, constructed the unpolarized (helicity-summed) cross-section, and then performed 
an integration over the orientation angles of the lepton pair (in the virtual-photon rest frame). 
After accounting for the different versions of dimensional regularization used [15], we have verified 
numerically that the two sets of results agree, for both the four-quark and two-quark-two-gluon 
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final states.''" Also, the squared matrix elements in refs. [13,14] have been incorporated into a NLO 
program for four-jet fractions and shape variables by Nagy and Trocsanyi [16], and their numerical 
results for the four-jet fractions confirm the corresponding results of ref. [12]. 

The amplitudes we present retain all correlations between the daughter or parent leptons of 
the vector boson, and the colored partons in the process. Such correlations are important for 
computations that take into account experimental constraints. For example, in the production 
of a W along with jets at a hadron collider, followed by the decay W — > £i>e, the longitudinal 
component of the W momentum cannot be observed, only that of the decay lepton, and the latter 
should be isolated from the hadronic jets in order for the event to pass detector cuts. In the case 
of jet production in deep-inelastic scattering, the orientation of the entire event with respect to 
the detector is dictated by the lepton scattering angle as well as the square of the virtual photon 
four-momentum. 

Recent years have seen a number of technical advances in the computation of one-loop am- 
plitudes. These advances have made possible the calculation of all one-loop five-parton processes 
[17,18,19], as well as a number of infinite sequences of loop amplitudes [20,21,22,23,24]. Many of the 
techniques used in the present calculation have been reviewed in ref. [25]. The rather complicated 
six-body kinematics encountered here necessitate further techniques for removing certain spurious 
singularities from the amplitudes. The removal of spurious singularities is essential in order to find 
(relatively) compact final expressions, and also plays a role in improving the numerical stability of 
the results. 

The general strategy employed in this paper is to obtain amplitudes directly from their an- 
alytic properties instead of from Feynman diagrams. In particular, we use the constraints of 
unitarity [26,22,23,27] and factorization [28,29,20,30], as summarized in ref. [25]. We construct the 
amplitudes by finding functions that have the correct poles and cuts in the various channels. The 
required pole and cut information is extracted from previously obtained amplitudes (tree ampli- 
tudes or lower-point loop amplitudes); manifest gauge invariance is therefore maintained at each 
step. This approach leads to compact expressions, especially when compared with those obtained 
from a traditional diagrammatic computation. In a Feynman diagram approach each diagram alone 
is not gauge invariant, and is often much more complicated than the sum over all diagrams. As a 
check, we performed a numerical evaluation of the Feynman diagrams at one kinematic point, and 
verified that their sum agrees with a numerical evaluation of our analytic results. 

The spinor helicity method [31] and color decompositions [32] are crucial to the success of 
this approach, because they simplify the analytic structures that must be computed. Because 

t We thank J.M. Campbell and E.W.N. Glover for assistance in comparing the results. 
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of the intricate analytic structure of the amplitudes it is rather non-trivial to remove spurious 
singularities that can appear in the amplitudes. (By a spurious singularity we mean a kinematic 
pole or singularity whose residue vanishes.) By evaluating the cuts appropriately we can prevent 
the worst of the spurious singularities from appearing. However, some of the spurious singularities 
are inherent in the amplitudes when they are expressed in terms of logarithms and dilogarithms. As 
we shall see, the spinor helicity method is quite useful for simplifying the terms containing spurious 
singularities. 

This paper is organized as follows. In section 2, we briefly review spinor helicity and color 
decompositions and provide formulas relating the full amplitudes to the 'primitive' amplitudes in 
terms of which the results are expressed. In section 3, we outline the construction of amplitudes 
from their analytic properties. Procedures for eliminating or simplifying spurious singularities 
are given in section 4. Sample calculations are given in section 5; in particular, examples are 
provided for cut constructions, rational function reconstructions, and simplifications via numerical 
analysis. The general structure of the primitive amplitudes including regularization issues is given 
in section 6. The results for the primitive amplitudes are collected in sections 7-12. In section 7 
the 'master functions', which are a set of functions which appear in multiple amplitudes, are 
given. Section 8 contains the amplitudes which are leading in the number of colors and flavors. 
Subleading-in-color contributions are contained in sections 9 and 10. Contributions with the vector 
boson coupled to a closed fermion loop are give in section 11; this includes both vector and axial- 
vector contributions. Simplified versions of the four quark amplitudes (which have been previously 
presented in ref. [10]) are given in section 12. Some concluding remarks in are given in section 13. 
There are a total of four appendices. Appendices I and II concern the evaluation of loop integrals 
and their associated spurious singularities. Appendix III lists some spinor-product identities that 
are useful for simplifying the spurious pole structure of kinematic coefficients. Appendix IV records 
the helicity amplitudes for e + e~ — ► qqg [7] in the same notation and conventions used in the paper; 
these amplitudes are useful because they appear in many collinear limits of the e + e~ — ► qqgg 
amplitudes. 



2. Brief Review of Basic Tools 

We shall present our results in terms of the spinor helicity method and SU(N C ) color decom- 
positions. The reader is referred to review articles [29] and references therein for details beyond 
our following brief review. 
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(2.2) 



2.1 Spinor Helicity 

We represent the gluon polarization vectors in terms of Weyl spinors | k ± ) [31], 

where k is the gluon momentum and q is an arbitrary null 'reference momentum' which drops out 
of final gauge-invariant amplitudes. The plus and minus labels on the polarization vectors refer 
to the gluon helicities. Our (crossing-symmetric) convention takes all particles to be outgoing, 
and labels the helicity and particle vs. antiparticle assignment accordingly. For incoming (negative 
energy) momenta the helicity and particle vs. antiparticle assignment are reversed. It is convenient 
to define the following spinor strings, 

(ij) = (K\kf) , M = (kt\kJ), 

W\3) = (K I h I k~) , m + m)\j) = (kj | (ft + f m ) | kj) , 

[i\lm ■ ■ ■ \j] = (kf | ftiftm • • • | kj) , 
(i\lm ■ ■ ■ \j) = (k~ | faftrn ■■■ | kf) , 

+ m){n + r) • • • \j] = {kf \ (ft + ^ m )(^ n + ft.) • ■ ■ | kj) , 
{i\(l + m){n + r) • • • \j) = (k~ | (ft + ^ m )(^ n + ft.) ■ ■ ■ | kf) , 

which is the notation we shall use to quote the results. In the last definition we take the | kf) to 
mean | kj) for an odd number of gamma-matrices in the string and | kf) when there are an even 
number. All the momenta ki are massless, kf = 0. Sometimes we will also use the notation 

(i\£ m \j)^(k-\£ m \kj), (2.3) 

etc., where t m is a loop momentum. The spinor inner products (ij), [ij] are antisymmetric and 
satisfy (ij) [ji] = 2ki ■ kj = Sij. In addition to s^- = (ki + kj) 2 we also define the three-particle 
invariants t^j = (ki + k 3 + ki) 2 . 

2.2 Color Decomposition 

It is convenient to decompose one-loop amplitudes in terms of group-theoretic factors (color 
structures) multiplied by kinematic functions called 'partial amplitudes' [32]. We present results 
for the general gauge group SU(N C ) (N c = 3 for QCD), and normalize the group generators in the 
fundamental representation, T a , so that Tr(T a T b ) = S ab . Color decompositions are obtained by 
rewriting the structure constants f abc as 

jabc = — 1 Tv([T a ,T b ]T c ) . (2.4) 
V2 



Then one applies the SU(N C ) Fierz identity 



(X 1 T a X 2 ) {Y 1 T a Y 2 ) = {X 1 Y 2 ) {Y 1 X 2 ) - ^(X 1 X 2 ) (Y l Y 2 ) , (2.5) 



where X i: Yi are strings of generator matrices T ai , in order to remove contracted color indices. 

Here we are interested in the amplitude Ae(lq, 2, 3, 4^; 5e, 6 e ), where legs 1,4 are the quark- 
anti-quark pair, legs 2, 3 are the gluon legs, and legs 5,6 are the lepton pair. We label the (outgoing) 
quark, anti-quark, electron and positron lines with subscripts q, q, e, and e, while the gluon lines 
do not have additional labels. The color decomposition of the tree-level contribution to As is 

X ree (l 9 ,2,3,4 9 -) = 2e 2 g 2 (-Q q + vl R vl R P z (s m )) £ (T»*WT°*w)> A^ e (l q , a(2), a(3), 4,) . 

(2.6) 

Here we have suppressed the 5, 6 labels of the electron pair, e is the QED coupling, g the QCD 
coupling, Q q is the charge of quark q in units of e, and the ratio of Z and photon propagators is 
given by 

*<"- .-Mi + «r,jft - (2J) 

where and are the mass and width of the Z. 

The left- and right-handed couplings of fermions to the Z boson are 



e _-l + 2sin 2 ^ e _2sm 1 w 

Vl ~ s\n20 w ' Vr ~ sm29 w ' 

q _ ±1 - 2Qg sin 2 9 W q _ 2Qi sin 2 6 W 

Vl ~ sin 26 w ' Vr ~ sin 26 w 



(2.8) 



where 9w is the Weinberg angle. The two signs in v q L R correspond to up (+) and down (— ) type 
quarks. The subscripts L and R refer to whether the particle to which the Z couples is left- or right- 
handed. That is, v q R is to be used for the configuration where the quark (leg 1) has plus helicity and 
the anti-quark (leg 4) has minus helicity, which we denote by the shorthand (1^,4^). Similarly, v q L 
corresponds to the configuration (1~,4^). Because the electron and positron are incoming in e + e~ 
annihilation, our outgoing-momenta notation reverses their helicities and particle vs. anti-particle 
assignment. Thus, v e R corresponds to the helicity configuration (5g ,6+) whereas w£ corresponds to 
the configuration (5g",6~). 



The one-loop color decomposition is given by 



^- loop (l„ 2, 3, 4,) = 2e 2 <? 4 { (-Q* + vl iR ^ R Vz{8 K )) 

x iV c ^(r-mr'O)^ A 6;1 (l q ,a(2),a(3),4«)+S a2a3 A 6;3 (l g , 4 g -; 2, 3) 

L <reS 2 

™f -1 



{T a 2T a 3) ^ + (r as T a 2) . u _ _ ^a, g.U ^ ;4 (1 ? , 4 ?S 2, 3) 



+ 



sin 29\y 



E ((^WT-'CS)).^ _ J-J^s A-(l„4 9 -;a(2),a(3)) 



<tgS 2 



+ _ J a 2 a 3 ^A-(l g ,4 ? -;2,3) 

(2.9) 

where Q l is the electric charge (in units of e) of the ith quark and rif is the number of light quark 
flavors. The partial amplitudes Aq. 4 , and A^ b represent the contributions from a photon or 
Z coupling to a fermion loop through a vector or axial-vector coupling. We take all quarks to be 
massless except the top quark. We assume that the top quark mass squared, m|, is larger than 
the other kinematic invariants in the process, and expand the fermion loop contributions in l/m 2 , 
keeping terms of order but dropping those of order 1/m 4 . In this approximation the top 

quark contribution to Aq. 4 vanishes (see section 11). On the other hand, in the axial vector channel 
isodoublet cancellations for massless quarks ensure that only the t, b isodoublet contributes to Af^ 
and Aq* 5 . There are also order l/mj contributions to A 6 -i. 

For convenience we also quote the color decomposition for the four-quark amplitudes. At tree 
level, we have 



-^6 ree (lgi 2q, 3q, 4, 



2eV 



(-<^ + vl iR vl R Pz(s 56 ))A^(l q ,2 Q ,3 Q ,4 q ) 
+ (~Q Q + vi,Rvi R V z (s 56 ))A% ee (3 Q ,4„ l q , 2 Q ) 



(2.10) 



while the one-loop decomposition [10] is 



^- loop (l„2 ,3 Q ,4,) 



2eV 



{~Q q + vl, R vl R Pz(s 56 )) [N c C A 6;1 (l q , 2 Q ,3 Q ,4 q ) + C C 4s ;2 (l„ 2g, 3 Q , 4,) 
+ + vl R vl R V z (s 56 )) [N c 6£ 6g A 6;1 (3 Q) 4„ 1„ 2q) + $ * A 6;2 (3 Q , 4„ 1„ 2q) 

+ Sdfc*<«»> ('* C " C)^;3(1 9 ,2 ,3 Q ,4 9 -) . 

(2.11) 

For the case of identical quark flavors (q = Q) see ref. [10]. We also include here contributions of 
order 1/m 2 from vacuum polarization loops to Aq-\ and Aq- 2 (which are very small at present e + e~ 
machines). The partial amplitudes A e -x and A e . 2 also appear in W + 2 jet production at hadron 
colliders; we leave the coupling constant and mass conversions as an exercise. 

The virtual part of the next-to-leading order correction to the parton-level cross-section is 
given by the sum over colors of the interference between the tree amplitude „4,g ree and the one-loop 
amplitude „4g _loop . Using the color decompositions (2.6) and (2.9), and the Fierz rules (2.5) the 
color-sum for e + e~ — ► qqgg in terms of partial amplitudes is, 

£ [^AInlo = 8eV (N 2 C - l)Re( (-Q« + v^ R v% R V z (s, & ))A^*{l q , 2, 3, 4 5 -) 



colors 



{~Q q + vI, r vIr V z ( S5 e)) [(N 2 C - l)A 6;1 (l q , 2, 3, 4,) 

-A 6;1 (l 9 ,3,2,4 9 -) + ^ 6;3 (l 9 ,4 9 -;2,3)] 

+ E("^ + H,R« + v R )Vz(s 56 )) [N c - — )^. 4 (1„4,-; 2, 3) 



i=l 



+ iSL p ^[(^-^)^ 4 ^' 3 ' 



sin 2^vf 

- Aa«(1„4,-;3,2) + ^^ 5 (1„4,-;2,3)] | + {2 - 3} , 

(2.12) 

The corresponding formula for e + e~ — ► qqQQ is straightforward to obtain, but lengthier, and so 
we omit it here. 

2.3 Primitive Amplitudes 

One may perform a further decomposition of the partial amplitudes in terms of gauge invariant 
primitive amplitudes [19,10]. Primitive amplitudes are gauge-invariant objects from which we can 
build partial amplitudes. They can be defined as the sum of all Feynman diagrams with a fixed 
cyclic ordering of the colored lines, with a definite routing of the fermion lines, and with vertices 
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that are given by color-ordered Feynman rules [29]. These differ from ordinary Feynman rules in 
that they have been stripped of the usual color factors. Here we choose not to fix the cyclic ordering 
of the colorless lepton pair with respect to the colored partons. In this case the lepton pair plays 
no role in the color structure, and so the equations expressing the partial amplitudes A & -i and ^6:3 
as sums of primitive amplitudes are identical to those derived in ref. [19] for one-loop two-quark 
two-gluon amplitudes — see eqs. (2.13) below. 



The main utility of primitive amplitudes as compared to partial amplitudes (i.e. the coefficients 
of a particular color structure) is that their analytic structures are generally simpler because fewer 
orderings of external legs appear. They also contain somewhat more color information, making it 
more straightforward to convert the amplitudes to those for other processes (for example, to replace 
gluons by photons, or quarks by gluinos). 



For the e + e~ — > qqgg amplitudes, 'parent diagrams' for each gauge invariant class are depicted 
in fig. 1. By a 'parent' diagram we mean a diagram from which all other diagrams in the set can 
be obtained via a continuous 'pinching' process, in which two lines attached to the loop are brought 
together to a four-point interaction — if such an interaction exists — or further pulled out from the 
loop, and left as the branches of a tree attached to the loop. The cyclic ordering of external legs is 
always preserved by pinching. Because we do not fix the ordering of the lepton pair with respect 
to the partons, primitive amplitudes with an external gluon cyclicly adjacent to the vector boson 
(j*,Z) have more than one parent diagram. For example, A e (l q , 2, 3, 4 q ) (fig. lb) has only one 
parent diagram, while ^(lq, 2, 3 q , 4) (fig. lc) has two parent diagrams, and A 6 (l q , 2 q , 3, 4) (fig. Id) 
has three. Note that these primitive amplitudes are different functions, not the same function with 
a different permutation of the arguments: A 6 (l q , 2 q , 3, 4) / Ae(l q , 2, 3, 4 q ). 
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(a) ^(l,, 2,3,4,- 



(h) A 6 (l g , 2,3,4,- 



c Ac 




(l„2,3,-,4) ^/\<^^ 



^(l„2 ? -,3,4) 



(e) ^ s,/ (l„2 ? -,3,4) 




(/) ^ x (l 9I 2,-,3,4) 



Figure 1. Parent diagrams for the various two quark and two gluon primitive amplitudes. 
Straight lines represent fermions, curly lines gluons, and wavy lines a vector boson (7* or Z). 

It turns out to be useful, in diagrams of the type shown in figs, lb, c and d, to replace a 
gluon loop contribution with two separate contributions, that of a scalar and that of the difference 
of a gluon and scalar, as shown in fig. 2. This separates the gluon loop contribution into two 
gauge-invariant pieces. As we shall discuss in section 3 this separation is advantageous because of 
differing analytic properties of the two pieces. Terms where a scalar replaces a gluon are labeled with 
a superscript 'sc', while those with the difference of gluon and scalar are labeled with a superscript 
'cc'. (As we shall discuss in section 3.1 the 'cc' terms are 'cut-constructible' meaning that they can 
be constructed solely from four-dimensional cuts.) Furthermore, terms proportional either to the 
number of scalars* n s (which vanishes in QCD), or fermions rif, are separately gauge invariant so we 



As in refs. [17,19], each scalar here contains four states (to match the four states of Dirac fermions) so that n s 
must be divided by two for comparisons to conventional normalizations of scalars. 
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also separate these out explicitly. (Due to Furry's theorem [33], or charge conjugation invariance, 
these terms appear only in Aq-i.) 






\ 

Figure 2. The contribution from a gluon in the loop is separated into the difference of a 
gluon and scalar, plus a scalar. 



Thus we take the decomposition of the partial amplitudes in terms of the primitive amplitudes 



to be 



4 6 ;i(l„ 2, 3, 4 q ) = A 6 (l q , 2, 3, 4 q ) - — A 6 (l q , 4 q , 3, 2) 



+ 



n. 



n f 



A s 6 (l q ,2, 3, 4 q ) - jt-4{l q , 2, 3, 4 q ) + ^A(l q , 2, 3, 4 q ) , 



(2.13) 



A 6;3 (l 9 , 4 9 -; 2, 3) = A 6 (l„ 2, 3, 4 g -) + A 6 (l„ 3, 2, 4 q ) + A>(l g , 2, 4 g -, 3) + A 6 (l q , 3, 4 g -, 2) 

+ A 6 (l q ,4 q ,2,3) + A 6 (l q ,4 q ,3,2) , 
Al. 4 (l q ,A q ; 2, 3) = -A v 6 s (l q ,4 q , 2, 3) - A v 6 f (l q ,4 q , 2, 3) , 
A-(l 9 ,4 q -;2,3) =^ x (l 9 ,4 9 -,2,3), 
^ 5 (l„4 5 -;2,3) = A^ s \l q ,4 q ,2,3). 

We have decomposed the fermion loop in Aq. 4 into a scalar loop piece Aq s and an additional piece 
A V J '. As mentioned above, we perform a further decomposition of the A 6 into cc and sc pieces, as 
depicted in fig. 2, 

A 6 (l„ 2, 3, 4 g -) = A-(l 9 , 2, 3, 4,-) + ^ c (l 9 , 2, 3, 4,-) , 

A 6 (l 9 , 2, 3,-, 4) = 2, 3,-, 4) + ^ c (l 9 , 2, 3 9 -, 4) , (2.14) 

A 6 (l q , 2 q , 3, 4) = AZ(l q ,2 q , 3, 4) + ^ c (l 9 , 2 q , 3, 4) . 

Finally, A\ gives the top quark contribution to Aq-i, through order 1/mf. 

We choose a set of helicity amplitudes from which all others may be obtained by applying the 
discrete symmetries of parity and charge conjugation. Parity reverses all external helicities in a 
partial amplitude; it is implemented by the "complex conjugation" operation, which substitutes 
00 [Ij], F° r the axial- vector fermion loop partial amplitudes one must multiply 

by an additional overall minus sign. Charge conjugation changes the identity of a fermion to an 
anti-fermion and vice-versa. These operations allow us to fix the helicity of the electron and the 
quark to be positive, and the positron and anti-quark to be negative. In addition, if the two gluons 
have the same helicity, we can fix that common helicity to be positive. 
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For the four-quark amplitudes, a similar use of charge conjugation and parity reduces the 
independent partial amplitudes to A 6;i (l+, 2~t, 3~, 4^~, 5J, 6+) and A 6;i (l+ , 2q, 3+, 4^ , 5jT, 6+). The 
formulas for the four-quark partial amplitudes, A 6;i (1^, 2q, 3^, 4^~), in terms of primitive amplitudes 
are [10] 



A 6;1 (l+ 2+,3- 4r) = A++(l,2,3,4) - A (^++(1, 2,3,4) + A+"(l, 3, 2, 4)) + ^^(2,3,1,4) 

+ ^^^' ++ (l,2,3,4)-^4' ++ (l,2,3,4) + i-4' ++ (l,2,3,4), 

A 6 ;2(l+,2+,3 Q ,4r)=^+-(l,3,2,4) + ^(^+-(l,3,2,4)+^++(l,2,3,4)) - ^<(2, 3, 1, 4) 

" ^p^<++(l,2,3,4) + ++(1,2, 3,4) - ^'++(1,2,3,4), 



A 6;3 (l+,2t,3-,4 g -) = ^ x (l,4,2,3), 

(2.15) 



and 



^6;i(l+,2^,3+ 4r) = A+-(l,2,3,4)- ^(^+"(1,2,3,4) +A++(1,3,2,4)) -^^'(3,2,1,4) 

+ ^^+"(1,2,3,4) - ^4 + -(l,2,3,4) + ^-4' + " (1,2, 3, 4), 

^ 6;2 (l+,2^,3j,4r) = ^++(l,3,2,4) + ^(A 6 ++ (l,3,2,4)+A 6 + -(l,2,3,4)) + ^<(3,2,1,4) 

_ V^L A ^ +" (i, 2 , 3 , 4) + ^4 + " (1, 2, 3, 4) - ±4- + " (1, 2, 3, 4) , 

4s ;3 (l+ 2g,3+, 4r) = -^-(1,4,3,2). 

(2.16) 

Although color decompositions do not depend on the helicity choices, the sign differences in these 
equations appear because we have used symmetries of the four-quark primitive amplitudes to reduce 
the number of independent ones required. In fig. 3 we display the parent diagrams associated with 
each of the primitive amplitudes 
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Figure 3. Parent diagrams for the four-quark primitive amplitudes. In each case the vector 
boson can appear on either side of a gluon line that is attached to the same quark line. 

3. Analytic Construction of Amplitudes 

In this section, we review the construction of one-loop amplitudes starting from their analytic 
properties. The two analytic properties we shall use are the determination of imaginary parts 
by Cutkosky rules, and factorization on particle poles. These properties of amplitudes have, of 
course, played an important role in field theory for many decades; the recent development which 
we focus on here is the ability to obtain, efficiently, complete amplitudes with no subtractions or 
ambiguities. These techniques, reviewed in ref. [25], have been applied to obtain results for both 
nonsupersymmetric and supersymmetric maximally helicity violating n-point amplitudes [20,22,23], 
and more recently for multi-loop N = 4 supersymmetric four-point amplitudes [24]. 

3.1 Cutting Rules 

Cutkosky rules [26,34] allow one to obtain the imaginary parts of one-loop amplitudes di- 
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rectly from products of tree amplitudes. (By imaginary parts we mean absorptive parts, that is 
discontinuities across branch cuts.) We apply Cutkosky rules to amplitudes instead of diagrams 
because amplitudes, being gauge invariant, are simpler. In the channel with momentum squared 
K 2 = (k mi + k mi+ i + • • • + k m2 ) 2 the cut of an amplitude is (see fig. 4) 



4T lo °P(l,2,...,n) 



K 2 cut 



= i J d D LIPS(-£!,£ 2 ) A tree (-£ 1 ,m 1 ,...,m 2 ,e 2 ) x A tTee (-£ 2 , m 2 + 1, . . . , m x - 1, h) , 

(3.1) 

where the integration is over D-dimensional Lorentz-invariant phase space with with momenta —l\ 
and l 2 for the intermediate states. For this channel K 2 is taken positive and all other kinematic 
invariants are taken negative. 




m 2 + 1 



- 1 

Figure 4. A cut amplitude, with momentum K 2 flowing across the cut. The lines represent 
gluons, scalars or fermions. 

The Cutkosky rules determine imaginary parts of the amplitudes. Dispersion relations are 
conventionally used to reconstruct real parts from imaginary parts. Here we bypass dispersion rela- 
tions; instead we replace phase-space integrals with cuts of unrestricted loop momentum integrals, 



4r loop (i,2,--- 



n 



I 



K 2 cut 

d D e 1 



(2n)° 



1 



^2 



A^i-h, mi, ...,m 2 ,e 2 )i2 A tree (-£ 2 ,m 2 + 1, . . . , m x - 1, h) -j 



K 2 cut 

(3.2) 

Whereas eq. (3.1) includes only imaginary parts, eq. (3.2) contains both real and imaginary parts. 
As indicated, eq. (3.2) is valid only for those terms with a if 2 -channel branch cut; terms without 
such a branch cut may not be correct. A useful property of this formula is that one may continue 
to use on-shell conditions for the cut intermediate legs inside the tree amplitudes without affecting 
the result. Only terms containing no cut in this channel would change. We are able to avoid the 
use of dispersion relations because we have additional information, namely that the reconstructed 
analytic functions are given by Feynman loop integrals. 

A similar equation holds for every branch cut. If one now combines all cuts into a single 
function having the correct cut in each channel, one obtains the full amplitude — up to the possible 
addition of a rational function, i.e. a function having no cuts at all, if one approximates the (4 — 2e)- 
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dimensional cuts by their four-dimensional limits (see discussion below). The full amplitude can 
also be expressed as a linear combination of various types of basic one-loop integral functions, 
multiplied by rational function coefficients. Many of the integral functions, for example scalar box 
integrals (which depend on the momenta of four external legs), have cuts in more than one channel. 
The coefficients of those integral functions as extracted from cuts in different channels must agree, 
and this provides a strong consistency check on the reconstructed amplitude. 

In general, it is convenient to take the tree amplitudes on either sides of the cuts to be four- 
dimensional. This is natural in the helicity formalism, which implicitly assumes that momenta are 
four-dimensional. On the other hand, we wish to regulate the ultraviolet and infrared divergences 
by letting D = 4 — 2e in the loop integral (3.2). When the (4 — 2e)-dimensional cut momenta 
are replaced by four-dimensional momenta an error may occur. Although the (— 2e)-dimensional 
parts are implicitly of 0(e), if the associated integral has an ultraviolet pole in e an 0(e°) rational 
function may remain. (As discussed in refs. [22,23], infrared poles never give rise to such rational 
contributions.) Despite this seeming error, a complete reconstruction of an amplitude is possible 
even when the 0(e) parts of the cut momenta are dropped, if the amplitude satisfies a certain 
power-counting criterion; we call such amplitudes cut-constructible [23,25]. The power-counting 
criterion is that the n-point integrals appearing in the amplitude should have (for n > 2) at most 
n — 2 powers of the loop momentum in the numerator of the integrand; two-point integrals should 
have at most one power of the loop momentum. (By an n-point integral we mean an integral 
with n propagator factors in the denominator, as in equation (1.1).) Cut-constructible amplitudes 
are composed of a restricted set of integral functions, and sufficient information exists from the 
four-dimensional cuts to determine the coefficients of each such function. These integral functions 
automatically include the cut-free rational functions in the amplitude [23] . 

The full e + e~ — ► qqgg and e + e~ — > qqQQ amplitudes are not cut-constructible. Diagrams 
containing a closed fermion, scalar or gluon loop can have up to n powers of the loop-momenta in the 
numerator of an n-point integral; other diagrams typically have up to n— 1 powers. However, one can 
split the amplitudes into 'scalar' contributions, plus additional terms which are cut-constructible; 
in fact, we have already performed such a decomposition in eqs. (2.13) and (2.14). In the case of a 
closed fermion loop, we wrote the fermion loop as the negative of a scalar loop, plus a second term 
(which is the contribution of an A = 1 supersymmetric chiral multiplet). As is apparent in a second- 
order formalism for fermions, as reviewed in ref. [25], the latter contribution is cut-constructible 
since the leading two powers of loop momentum cancel. 

In the case where an external quark line attaches directly to the loop, for example the leftmost 
diagram in fig. 2, an n-point integral has a maximum of n — 1 powers of the loop momentum. If 
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one subtracts from this diagram an identical diagram but with the gluon replaced by a scalar (also 
shown in fig. 2), suitably adjusts the scalar-fermion Yukawa coupling, and works in background- 
field Feynman gauge [35], then the leading loop-momentum terms cancel. Thus the difference is 
cut-constructible [23] . 

In both cases, although the scalar contributions are not cut-constructible, they are simpler in 
some respects than the full amplitudes; for example, certain cuts vanish in the scalar contribution, 
but not in the full amplitude. They also have spurious singularities of different degree from the 
cut-constructible terms. To determine the rational functions for such amplitudes, which do not 
satisfy the power-counting criterion (but contain only massless particles), one may compute to one 
higher power in the dimensional regularization parameter e [36,27,25]. However, a more convenient 
approach here is to ignore all 0(e) contributions, and instead use the amplitudes' factorization 
properties to reconstruct their rational functions. 

3.2 Factorization 

The properties of tree-level QCD amplitudes as kinematic invariants vanish have been pre- 
sented in various reviews [29]. The corresponding one-loop factorization properties have also been 
extensively discussed [20,30,25], so here we will briefly review only the salient features necessary 
for obtaining the rational function parts of amplitudes. 

For amplitudes with six- and higher-point kinematics, the properties of amplitudes under 
factorization when a kinematic invariant vanishes are in general sufficiently powerful to obtain the 
rational function parts of amplitudes. Although factorization is complicated by the appearance of 
infrared divergences, nevertheless, as any kinematic variable vanishes, one-loop amplitudes have 
a universal behavior quite similar to that of tree-level amplitudes [22,30]. If one finds a function 
which obeys the proper factorization equations in all channels, one has an ansatz for the rational 
function part of an amplitude. Although no proof of the uniqueness of such a construction has 
been presented, for six- and higher-point amplitudes no counterexample is known. (For a five- 
point amplitude counterexample see ref. [37].) Factorization can be a particularly efficient way to 
obtain the rational function terms; it avoids the need to perform loop integrals. This complements 
the efficiency of the cut-construction technique for obtaining the logarithms and dilogarithms. 
(Resorting to Feynman diagrams to obtain analytic expressions for the rational function parts is not 
satisfactory because such pieces tend to have the most complicated diagrammatic representation: 
they are associated with the maximal power of loop momenta and the largest number of diagrams.) 

Of particular utility are the two-particle factorization properties, depicted in fig. 5. As two 
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momenta become collinear the amplitudes behave as [20,22 
A i-ioo P «ll> £ f Split^V", b x ») A£T{. ■ ■ (a + b)- 

A=± ^ 



(3.3) 

+ Split L- loop (a\ b x ») AT.\(. ■ ■ (a + b) x . . .)) , 



where k a — > zP and A;;, — ► (1 — z)P, with P = k a +kb, P 2 = s a (, — ► 0. The helicity of the intermediate 
parton P is labeled by A. The tree and loop splitting amplitudes, Split and Split l_^ loop , behave 
as l/y/Sab in this limit. A complete tabulation of the splitting amplitudes appearing in one-loop 
computations in massless QCD has been given in refs. [22,19,30]. Given the n — 1 point amplitude 
and splitting amplitudes (or 'factorization functions' in multi-particle channels), eq. (3.3) provides 
an extremely stringent check since one must obtain the correct limits in all channels. A sign or 
labeling error, for example, will invariably be detected in some limits. 






+ n 



a ■ 

Figure 5. A schematic representation of the behavior of one-loop amplitudes as the momenta 
of two legs become collinear. 

The physical poles that can appear in any massless amplitudes are square-root singularities in 
two-particle channels, or single poles in multi-particle channels [29,25], 

11 11 1 . a , 

JT^~—, 77^-—, — • (3.4) 

Any other kinematic pole-type singularity that appears in individual terms of an amplitude must 
be spurious; that is, the residue of the pole must vanish for the full amplitude. 

In this paper we use factorization to construct ansatze for the rational function parts of the 
amplitudes. Although the construction of such an ansatz involves a certain amount of guesswork, 
the procedure can be systematized somewhat. In general, the rational functions contain non- 
removable spurious singularities (see section 4). In the full amplitude, this singular behavior cancels 
against singular behavior in terms that have logarithms and dilogarithms, and so it can be inferred 
from the information provided by the cuts. Therefore one can readily introduce rational function 
terms that account for all, or at least most, of the spurious singularities. Subtracting these terms 
from the full rational function ansatz leads to an ansatz for the remainder which is (largely) free of 
spurious singularities. At this stage it is relatively straightforward to proceed channel by channel, 
adding terms to the ansatz that correctly reproduce the desired singular behavior in each channel, 
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using the factorization information provided by lower-point amplitudes. The last few channels 
are simplest, because by then the remaining terms have very few singularities left. For six-point 
kinematics this procedure invariably gives the correct result, as we have verified by numerical 
comparisons to Feynman diagram computations. We determine the rational function terms for a 
simple example in section 5.2. 



4. Spurious Singularities. 

In this section we discuss the procedure for removing or at least greatly simplifying spurious 
kinematic singularities appearing the amplitudes. This is essential in order to obtain compact 
results. The presence of large numbers of such singularities leads to unwieldy results which tend to 
be unsuitable for use in jet programs since they are numerically unstable. Although the expressions 
encountered in the calculation of the cuts via eq. (3.2) are rather compact when compared to a 
direct Feynman diagram calculation, the amplitudes are sufficiently intricate that even a small 
number of spurious singularities can seriously hinder attempts to obtain compact results. 

Spurious singularities fall into two categories, removable and non-removable: After expressing 
the amplitudes in terms of logarithms and dilogarithms there are spurious singularities that can be 
removed from the amplitudes by algebraic simplification, and ones that cannot be removed. As a 
trivial example, consider the functions 

1-x 2 Inx lnx + l-x 

/l = T^7 ; /2 = T^' h = {i- x y ' (41) 

where the poles at x = 1 have vanishing residues. The first function can be re-expressed as 
fi = (1+x) so the spurious pole at x = 1 is removable. However, the spurious poles in / 2 and fs are 
not removable if we wish to express the function in terms of a logarithm. One may, of course, define a 
new set of functions, of which / 2 and f-j are examples, which absorb the spurious poles. Indeed, this 
is the role of the Lj and Ls^ functions [17] defined in appendix II. (Generalizations of such functions 
have been presented in ref. [38].) In the functions x is a ratio of kinematic invariants; for 
example, x = S23A123 arises from the Gram determinant Ag 2,5 - 1 in eq. (1.12). In Feynman diagram 
calculations both removable and non-removable singularities occur. The removable singularities 
are an artifact of the integral reduction techniques employed, but the non-removable ones are an 
inherent part of the amplitudes when they are expressed in terms of logarithms and dilogarithms. 

The complications arising from spurious denominators follow largely from dimensional analysis 
and the fact that they carry positive dimensions. Their appearance implies that the numerators 
must have compensating powers of momenta. In a six-point amplitude there are five independent 
momenta, leading to a substantial proliferation in the number of possible numerators. As an 
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example, consider a five-point tensor integral with four powers of loop momentum in the numerator 
encountered in a Feynman diagram evaluation of e + e~ — > qqgg. If this integral were evaluated 
by conventional means using a Passarino-Veltman reduction [39], summarized in appendix 1.1, one 
would have up to four inverse powers of the pentagon Gram determinant A 5 defined in eq. (1. 10). 
As we shall show below, these spurious singularities can always be removed from the amplitudes 
of this paper, but if they appear in intermediate expressions, their removal is an arduous task. 
The appearance of a spurious Ag 4 denominator implies that the numerators must contain an 
appropriate polynomial to cancel the poles. Since the terms in the numerator of a brute force 
calculation appear in a seemingly haphazard pattern, one must deal with the order of 22 4 ~ 10 5 
terms to remove the spurious singularities in A 5 . Moreover, the various spurious singularities can 
get tangled together in rather intricate ways. Our goal is to eliminate those spurious poles which 
do not belong in the amplitudes and to simplify those which are are inherently associated with the 
logarithms and dilogarithms. 

4.1 Types of Spurious Singularities Appearing in the Amplitudes 

The simplest type of spurious singularities are unphysical poles in the kinematic variables s^- 
and tijk- Even at tree level the kinematic singularities in the amplitude can be non-trivial. For 
example, consider the tree amplitude Ag ree (l+, 2 + , 3~, 4^~) given in eqs. (8.8) and (8.9). In the 
first form (8.8) the amplitude exhibits the proper poles in the £123 and £234 three-particle channels, 
but the behavior as S23 or s 56 vanish is not manifest — the apparent full poles in these variables 
actually cancel down to square-root singularities. The second form (8.9) of the amplitude exhibits 
the proper square-root singularities in all two-particle channels, at the expense of more obscure 
behavior in three-particle channels. 

At one loop, the situation is greatly complicated by the large variety of spurious singularities 
that arise from loop integrals. The parent diagrams for the Aq primitive amplitudes in figs. 1 and 3 
require the evaluation of pentagon integrals where all internal lines are massless, and all external 
legs are massless except for the leg composed of the lepton pair 5-6, which has invariant mass 
S56- Besides this one-mass pentagon integral, there are a number of different types of box integrals, 
where either one or two external legs are massive. These integrals may appear directly in a diagram 
or a term in a cut evaluation, or they may be generated in the reduction of pentagon integrals, 
as reviewed in appendix I. Similarly, we find triangle integrals with one, two or three external 
massive legs, as well as bubble (two-point) integrals. All these integrals have associated with them 
different kinematic factors, which can appear in the denominators of coefficients of logarithms and 
dilogarithms, and whose vanishings correspond to separate spurious singularities. 

In appendix I we summarize some of the standard integral reduction methods [39,40,41,42], 
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and their associated spurious singularities. As discussed in the appendix, when evaluating tensor 
n-point loop integrals one obtains denominators containing (minus) the Gram determinants 

A n = - det(2Ki ■ Kj) , i,j = 1,2,... ,n-l, (4.2) 

where the are external momenta or sums of external momenta. Other kinematic denominators 
which appear are the determinants 

det(%), i,j = 1,2,..., n, (4.3) 

where the symmetric matrix Sij is given by 

S u = 0, S ij = -±(K i + K i+1 + --- + K j _ 1 ) 2 , iovi<j. (4.4) 

In appendix I we collect the explicit forms of the various determinants that can appear as poles 
in the amplitudes. Particularly important are poles in the three-mass triangle Gram determinant 

A 3 = A 3 2 ' 4) = s\ 2 + s| 4 + Sgg - 2 S12 s 34 - 2 s 34 s 56 - 2 s 56 S12 , (4.5) 

in the spinor strings 

(1|(2 + 3)|4), (4|(2 + 3)|1), (4.6) 

and in objects related to these by permutations of the external legs 1, 2, 3, 4. The spinor strings (4.6) 
vanish in the kinematic configuration where ^2+^3 = aki+bk^ where a and b are arbitrary constants. 
We call this configuration 'back-to-back' because in the center-of-mass frame for particles 1 and 
4, viewed as incoming, the outgoing three-momenta ki + ks and A; 5 + fc 6 must be parallel to k\ 
and k 4 . In many cases, these singularities cancel only after taking into account the behavior of 
the dilogarithms and logarithms that appear, in analogy to the behavior of f 2 and / 3 in eq. (4.1). 
The appearance of the three-external-mass triangle Gram determinants (4.5) and the back-to-back 
singularities (4.6) explains to a large extent the significant increase in complexity of the amplitudes 
presented in this paper, as compared to the massless five-parton amplitudes [17,18,19]. It is essential 
to simplify terms containing these singularities if we wish to obtain (relatively) compact expressions. 

4.2 Integral Reductions 

We now describe techniques we used to help minimize the algebraic complexity of intermediate 
expressions when evaluating a cut amplitude (3.2). In particular, these techniques prevent the ap- 
pearance of the pentagon Gram determinants, which are by far the most noxious of the unwanted 
determinantal denominators. The same techniques apply just as well to Feynman diagrams, al- 
though in evaluating a cut one generally begins with a much more compact expression, making it 
simpler to keep its size small. 
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An important aspect to the calculations performed in this paper is the use of a helicity basis 
for both quark and gluon external states. This basis simplifies general gauge theory amplitudes, 
as reviewed in refs. [29,25]. (To make effective use of the helicity formalism at the loop level it 
is important to use a compatible regularization scheme, such as the FDH scheme [43], which at 
one-loop has been shown to be equivalent to a helicity form of dimensional reduction [44,15].) The 
spinor helicity method also leads to a useful procedure for evaluating tensor integrals. The basic 
observation is that certain combinations of Lorentz invariant products, such as S56S23 — £123*234, 
which are destined to appear in the denominators of amplitudes for reasons discussed in the previous 
subsection, and which cannot be factored in terms of Lorentz invariants, can be factored in terms 
of spinor strings such as (1|(2 + 3)|4), as shown in eq. (1.16). By multiplying and dividing by such 
spinor square roots, and by further manipulating the spinor strings in the numerator of the loop 
momentum integral for a cut, one can extract inverse propagators. These factors cancel propagators 
in the denominator, leaving behind much simpler lower-point integrals to evaluate. This general 
strategy is similar to the more conventional Passarino-Veltman reduction in that it expresses tensor 
integrals as linear combinations of simpler integrals, but it differs in its economical use of expressions 
which already appear in the loop momentum integrands. The strategy has already been applied 
to a number of cases [22,23,18,27,45,24]. The amplitudes of this paper have rather complicated 
kinematics, and so it is very important to exploit factorization relations such as (1.16), in order to 
appropriately arrange the spinor strings and thus maximize the simplifications. 

When evaluating a cut (3.2) one typically encounters a pentagon loop momentum integral of 
the form 

d^p {a\ti\b){c\lj\dj--- 



I 



j ' ) — , I 2, rl rl rl rl rl (4-7) 
V / 1 2 3 4 5 

where = p — k\ — ■ ■ ■ — is the momentum of the ith loop propagator. The kinematic 
configuration of this pentagon integral is shown in fig. 6. By multiplying and dividing by (b\d\c), 
we may convert the numerator factors into a single spinor string containing both i i and jtj, 

(a\l t \b)WM = WM^\b) . (4.8) 

We may then extract inverse propagators by commuting l i and jtj towards each other. In this way 
we generate terms of the form 

2 kj • £ m , 2 £i • £j . (4*9) 

For this to be a useful rearrangement, these dot products must be expressible as sums of inverse 
propagators and external kinematic variables. This requirement dictates some care in deciding 
which pairs of strings (a\£i\b), etc., to work with, and which momentum d to use in the string (b\d\c) 
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that one multiplies and divides by. An inappropriate choice can lead to a large algebraic expression 
without any reduction in powers of loop momenta in the numerators of the integrals. 




Figure 6. The momenta associated with the pentagon integral appearing in the example. 
For example, consider a numerator containing the product 



(5|4|1)(4|4|6) 



(4.10) 



which occurs in the evaluation of the cut in the s^q channel for the amplitude vlg c (1+ , 2~ , 3 + , 4^ , 5J , 6+ ) . 
In this channel we may use the on-shell conditions l\ = and i\ = 0. Before explaining a good 
choice for forming a spinor string we mention first some choices which are not very helpful. For 
example, we might combine the spinors via 

1 ... .... , 1 



(5|4|1)(4|£ 5 |6) 



(5|^|1)(16) <6 + |4|4 + > 



<5|4 164>|4> 



(4.11) 



(16) x 1 11 101 ' (16) 

In this case, if we commute i x or /s P as t ^6 we obtain a term containing k 6 • l\ or fc 6 • £ 5 , neither 
of which can be expressed in terms of inverse propagators, since only fc 5 + ka appears in the 
loop propagators. Other choices, such as multiplying and dividing by (1|3|4), are better but still 
introduce unwanted spurious singularities in the amplitudes, which must be removed at later stages 
in the calculation. 

A much better choice is to multiply and divide by (1|(2 + 3)|4) = — (1|(5 + 6)|4), so that 

1 



(5|£ 1 |1)(4|£ 5 |6) 



■(51^1(5 + 6)4416), 



(4.12) 



(1|(2 + 3)|4) 

which will ensure that we commute i 1 and / 5 only with neighboring momenta. This choice is 
motivated by the appearance of this kinematic singularity in the scalar pentagon integral reduction 
formula (1.7) (with n = 5), after expressing the determinantal denominator (1.13) in the factored 
form (1.16). 

Once a numerator term is in a form where at least two t i are in the same inner product we can 
commute these terms toward each other. In particular, for the spinor string appearing in eq. (4.12) 
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we have 

(5|41(5 + 6)44|6) = ^(5|^il5|6) + ll(5\6£ 5 4\6) + (5|1^(5 + 6)£ 5 4\6) 

(4.13) 

= t\ [65] (51) (51^11) -^(56) [64]<4|4|6>, 

where we used 

<5|lfi(5 + 6)f 5 4|6) = (511^(4-^)4416) =0, (4.14) 

and 

2£ 5 -k 4 = (4 + k 4 ) 2 =£l, 2e 1 -k 1 = -4, (4.15) 

which follow from the s 56 -cut conditions l\ = l\ = 0. Since both terms in eq. (4.13) contain inverse 
propagators we have succeeded in reducing the pentagon integral to a sum of two box integrals. 
The rather clean simplification is due to our choice of multiplying and dividing by (1|(2 + 3)|4). 
Of course, not all cases are reduced as easily, but this example does illustrate the importance of 
choosing appropriate factors to multiply and divide by. 

Each inverse propagator appearing in a numerator cancels a propagator, leaving a lower-point 
integral with one less power of loop momentum; if there are any terms without an inverse propagator 
then they are also down by one power of loop momentum. Thus, by combining spinor strings and 
commuting pairs of / i toward each other, we can always reduce a pentagon integral with m > 1 
powers of loop momentum to a linear combination of pentagon and box integrals with at most m— 1 
powers of loop momentum. When one or no powers of loop momentum are obtained we may use 
the reduction formulas (1.7) and (1.9) (with n = 5), which are free of pentagon Gram determinants. 
In this way, we avoid encountering any A 5 s in the calculation. 



Integral 


Back-to-Back Singularities 


t r(3) f(3) 
J 5 ! A 4 i 1 4 


(1|(2 + 3)|4),(4|(2 + 3)|1) 


f r( 2 ) f( 2 ) f(2) 
1 5 > J 4 j 1 4 j 1 4 


(4|(1 + 2)|3),(3|(1 + 2)|4) 


f fW f(i) 


(3|(4 + 1)|2),(2|(4 + 1)|3) 


r(4) f(4) 


(1|(3 + 4)|2),(2|(3 + 4)|1) 



Table 1: The back-to-back singularities associated with reductions of each type of integral 
function. 

The best quantities to multiply and divide by when forming a spinor string are usually the 
back-to-back spinor products which would occur in the Passarino-Veltman integral reductions. 
In many cases, these singularities are not removable and appear in our final expressions. Using 
the results summarized in appendix I, we have collected in table 1 the back-to-back singularities 
associated with the reduction of each type of integral. This table provides guidance in actual 
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calculations as to which spinor products to introduce, although in some cases simpler alternatives 
are available (see below). Since the (7*, Z) is not colored, for a given color ordering it may appear 
with either cyclic ordering with respect to a cyclicly adjacent external gluon. This means that when 
the first four legs are ordered 1234 by the color flow, one may still have integral functions with three 
possible orderings for the six legs: 123456, 123564 and 125634. We denote the pentagon integrals 
corresponding to these three orderings by 1$, 1 5, and I 5 . The box integrals are similarly denoted 
by 1^ and 1^ where the label (i) indicates that the box is obtained from the pentagon 
by removing the propagator prior (in the clockwise ordering of legs) to the ith external leg. For 
convenience we have collected the integrals appearing in Table 1 in fig. 7. 




Figure 7. The integrals that can lead to spurious back-to-back singularities. 

Sometimes the pentagon integral 1$ appearing in a cut can be reduced to boxes without 
introducing the 'back-to-back' factor (1|(2 + 3)|4) or its complex conjugate (4|(2 + 3)|1). In such 
cases, the required factor one should multiply and divide by is either (23), or else its complex 
conjugate [23]. For example, the cut in the s$q channel for the amplitude Ag c (l+, 2^", 3~, 4 + , 5^ , 6+) 
contains a term with the factor 



(4+K 3 |2+>(3|4|4> = ^[4|^ 3 |4], 
25 



(4.16) 



which is easily reduced to 



jL M - 4) [43] (3|4|4) + (ll - I 2 ,) [42] (2|4|4) + £l [4|23|4]" . 



(4.17) 



After having reduced pentagon integrals to box integrals, the next step is to reduce the box 
integrals. Boxes with two adjacent massive legs, depicted in columns two and four of fig. 7, have 
the same kind of back-to-back singularities associated with them as does the pentagon integral 
(see table 1); therefore they can usually be reduced by multiplying numerator and denominator by 
the string (i\(k + l)\j) (or its complex conjugate), along the lines of eq. (4.12). Table 1 and fig. 7 
show that for this box reduction i and j should be the two adjacent massless legs. In contrast, 
for the other two types of box integrals that appear, boxes with only one massive leg and boxes 
with two diagonally opposite massive legs, the appropriate factor to multiply and divide by turns 
out to be a single spinor product, (ij) (or its complex conjugate), where now i and j represent the 

(5) 

two diagonally opposite massless legs. For a box integral with one massive leg, such as I\ , this 
is suggested by the factor of s\z = — (13) [13] in the Gram determinant in eq. (1. 11). For a 
box integral with diagonally opposite massive legs, such as I\ , the corresponding factor in A4 
is S14 = — (14) [14]. Even though the spinor string (1|(2 + 3)|4) also appears as a factor in it 
is never required for the reduction of this box integral. 

The appropriateness of these factors for simplifying numerators is due to the fact that they 
terminate with spinors corresponding to massless legs of the relevant integrals, not just the ampli- 
tude. The key to the previous pentagon reductions was having only massless legs of the pentagon 
integral interposed between two loop momenta in a single spinor string. At the level of box and 

lower-point integrals, fewer of the external legs of the integral are massless, because they can in- 

(2) 

stead be sums of massless legs of the amplitude. For example, in the box integral I\ in fig. 7, 
legs 3 and 4 are massless legs of the integral, but legs 1 and 2 only appear as constituents of a 
massive leg. If we had a string of the form ■■■1£ 1 ---, and we tried to commute l x past $1 , we would 
generate 1l\ -k\ = e\ — 1\. But is not a propagator for lf \ hence the commutation procedure 
fails to reduce the integral. (This is the same problem encountered when the momentum k& was 
introduced in the pentagon numerator (4.11).) 

In some cases, not enough of the massless spinors terminating the spinor strings in the numer- 
ator of an integral correspond to massless legs of the integral. For example, if the string 



( a \e 3 \b)(3\e 3 \c) 



(4.18) 



appears in the numerator of the adjacent two-mass box integral l\ , and none of a, b and c is equal 
to 3 or 4, then we are seemingly blocked from using the above procedures. One way to handle this 
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situation is to use the Schouten identities, 



(ij) (kl) = (il) (kj) + (ik) (jl) , 



(4.19) 



[ij] [kl] = [il] [kj] + [ik] [jl] 



to put more 'useful' momenta at the ends of strings. In the present case, we multiply and divide 
by [34] . Then we use 



and the last term can now be reduced further, after multiplying and dividing by (4|(1 + 2)|3). 

The trick of multiplying and dividing by spinor product factors stops working when one has 
too few massless legs, which here basically happens at the level of triangle integrals. Fortunately, 
the triangle integrals with one and two external masses do not generate terribly complicated ex- 
pressions even when reduced via a general ('brute force') formula, for example using Feynman 
parametrization. On the other hand, three-mass triangle integrals with two or three powers of the 
loop momenta inserted can generate quite lengthy formulae. Indeed such terms — coefficients of 
-^| m ( s i2) S34) S56), ln( 5~ ); etc. — account for much of the length of our final expressions. Part of 
the problem is that the three-mass triangle Gram determinant Ag 2 ' 4 ^ in eq. (1.12) — which 'belongs' 
in the various coefficients in some form — has mass dimension 4, yet apparently cannot be factored 
at all, even employing spinor strings. 

The cuts for one-loop six-point amplitudes can be divided into and tijk cuts, according 
to whether the momentum flowing across the cut is the sum of two or three external momenta. 
In general, the t^k cuts are much simpler to evaluate, largely because the three-mass triangle 
cannot appear in such a cut — it has cuts only in three Sij channels. One way to handle the 
more intricate Sjj cuts is to first evaluate a 'triple cut', where three loop propagators are required 
to be on-shell. For example, the S12-S34-S56 triple cut, depicted in fig. 8, can be defined by the 
conditions i\ = i\ = t\ = 0. Such triple cuts pick out those integral functions containing cuts in all 
three channels, in particular the three-mass triangle functions. The utility of considering such cuts 
is that instead of having to evaluate the phase-space integral of a six-point tree amplitude with a 
four-point tree amplitude (as one would for an s^j cut), one gets an expression where the six-point 
amplitude is replaced by the product of the two four-point amplitudes that it factorizes on. The 
full S56 cut (say) can then be written as a sum of the triple cut and a residual term which has 
no 1/^3 propagator, and therefore no three-mass triangle integral; that is, the three-mass triangle 
complications can be confined to the triple cuts. 



<a|£ 3 |&> [34] = <a|€ 3 |3> [64] + (a\£ 3 \4) [36] . 



(4.20) 



The string becomes 




£l) [64] (a\e 3 \c) - i\ [64] (a\3\c) + [36] (a|4|4)(3|4|c) , 



(4.21) 
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Figure 8. The kinematics of the triple cut. The cut lines are all on-shell. 

4.3 Numerical simplifications 

Even after employing the above spinor-product manipulations (among others) to help reduce 
the size of the expressions for cuts, one may still generate in intermediate steps complicated an- 
alytical expressions for the coefficients of If m (si2, S34, s$g), ln(— S12), and so on. It is possible to 
use numerical techniques to simplify such an expression, which is some rational function of the 
spinor products, if one has enough information about the analytic behavior of the coefficient. (The 
analytic information does not have to be manifest in the complicated expression.) The basic idea 
is to write down an ansatz for the complicated expression, as a linear combination of all possible 
kinematic terms that have the proper analytic behavior, where each term is multiplied by an as-yet 
unknown numerical coefficient. The more one knows about the analytic behavior of the coefficient, 
the fewer factors one has to put in the denominator of the ansatz, and (by dimensional analysis 
and combinatorics) the fewer the possible terms. Then one evaluates both the complicated expres- 
sion and the ansatz at a number of random kinematic points, which should exceed the number of 
unknown coefficients. This gives an over-determined set of linear numerical equations which can be 
solved for the unknown coefficients, which are required to be simple rational numbers. The solution 
can be checked by numerical evaluation at further kinematic points. An example of this procedure 
is provided in section 5.3. 

In practice we have been able to carry out this 'numerical simplification' procedure efficiently 
once the number of linearly independent terms in the ansatz for a given coefficient is reduced to 
about a hundred. To get a manageable number of terms like that, one generally needs to know 
about not just the physical factorization limits discussed in section 3.2, but also the leading spurious 
singularities of the coefficient. However, this latter information can be typically be inferred from 
simpler cuts which have already been performed. 
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For example, the coefficients of both I 3 m (si 2 , s 34 , s 56 ) and ln(— s 12 ) receive contributions from 
the adjacent two-mass box integral I\ , and thus they will typically have denominators containing 
(3|(1 + 2)|4). In the full amplitude the singularities in these terms as (3|(1 + 2)|4) — > cancel against 
singularities in the terms containing Ls 2 ^ h (s34, ii 23 ; s i2, s§§). The precise way this cancellation 
happens is known from the structure of the tensor box integrals [38]. In the present example we 
have, in the limit as (3|(1 + 2)|4) -> 0, 
[coefficient of /| m (si 2 , s 34 , s m )] 

r ffi • , f T 2mfc, , \i v / n <3|(1 + 2)|4)(4|(1 + 2)|3)(t 123 - t 12A )s l2 s 56 
-> [coefficient of Ls_ x (s 34 , t i23 ; s i2 , s 56 )J x (-1) — r , 

t 123 A 3 

[coefficient of ln(— s 12 )] 

r ffi • t f T 2mh( . n <3|(1 + 2)|4)(4|(1 + 2)|3)(t 123 - ^figjig 
-»■ [coefficient of Ls_ 1 (s 34 , tm; S12, s 5 6) I x (-1) — r . 

2s 34 t 123 Z\ 3 

(4.22) 

On the other hand, the coefficient of Ls^_™ /l (s 34 , t 123 ; s 12 , s 56 ) can be determined from the cut in the 
ti23 channel, which as we have mentioned is much simpler, typically only one or two terms in length. 
The known analytic behavior of the coefficients of l| m (si 2 , s 34 , s$q) and ln(— si 2 ) as (3|(1+2)|4) — > 
can also be verified numerically by choosing kinematics close to this 'back-to-back' singularity. 

Another denominator factor usually present in the coefficient of ln(— si 2 ) is (si 2 — ti 23 ). Again 
these terms can be inferred from the ti 23 -cut; this time the relation is through the simpler functions 
L (^) and Li(^u-). 

v — *123 ' v — *123 ' 

The singularity of ln(— si 2 ) as A 3 (s 42 , s 34 , S56) — > is related to that of /| m (si 2 , s 34 , s^e), but it 
cannot be related to a t^-fc cut, and has to be extracted from the leading loop-momentum behavior 
of the si 2 -s 34 -S56 triple cut. The numerical study of the A 3 — > limit is also subtler: In this limit 
the three vectors k\ + k 2 , k 3 + k 4 and k$ + ke all become proportional to each other; thus there 
is a simultaneous 'back-to-back' vanishing of (1|(3 + 4)|2), (3|(1 + 2)|4) and (5|(1 + 2)|6). Terms 
with A 3 in the denominator often have such vanishing factors in the numerator (or other more 
complicated ones — see eq. (1.17)), which mask the presence of I/A3. On the bright side, these 
numerator factors imply that individual kinematic coefficients are less singular as A 3 — > than a 
simple counting of denominator A 3 s would suggest, improving their numerical stability near the 
pole. 



5. Sample Calculations 



5.1 Evaluation of Cuts 

Here we describe the initial stages of evaluation of two different cuts, in order to illustrate 
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some of the features that are involved. 

The first cut we consider is the cut in the £234 channel of the amplitude Ag c (l+, 2^, 3~ , 4 + , 5 J , 6+). 
This amplitude has a scalar replacing the gluon in the loop, and so the cut propagator £ 2 is that of 
a scalar. (The configuration of external legs and loop momenta is shown in fig. 6.) The emission of 
the scalar from the quark line is via a Yukawa coupling that reverses the helicity of the quark line. 
Thus the required product of tree amplitudes, to be integrated over the two-body phase space for 
h + (-4), is 

P234 = 4 ree (l^(^),,(-4)^5 e -,6 e + )x4 ree ((-£ 2 ) s ,2r,3-,4+,(4)-). (5.1) 

These two five-point tree amplitudes are easily evaluated, up to overall signs (which can always be 
fixed at the end of the calculation, using e.g. a factorization limit), 

p _^ <^5) 2 „ [24] [l 2 A] 2 



(U2) (£ 2 £ 5 ) (56) [23] [34] [£ 2 £ 5 ] [£ 2 2] 
[24] (5|l 2 |4)(5|l 2 |l)(2|l 2 |4) 
[23] [34] t 234 (56) t\£\ 

In the second step we used the on-shell conditions i\ = i\ = to replace 1/ (1^ 2 ) with [£ 2 1] / (2£ 2 ■ ki) 
[£ 2 1] /£\-, and a similar replacement for 1/ [£ 2 2]. Notice that the propagator £\ is 'missing'. Hence 
no pentagon reduction is necessary in this example; the integral is already the box integral 1^ 
with two adjacent masses. (The missing propagator can be attributed to a supersymmetry Ward 
identity (SWI) [46]: The limit ^4 — > is also the collinear limit £5 \\ k^ for the tree amplitude 
A t 5 ree ((-£ 2 ) s ,2^,3-,4+,(4)-), in which it factorizes onto A\ vee ({-£ 2 ) s ,2g ,3~ , P~), which van- 
ishes by a SWI.) 

Table 1 and the structure of the numerator in eq. (5.2) suggest that we multiply and divide 
this expression by (1|(3 + 4)|2), and then commute the pair of / 2 s toward each other in the following 
string: 

(5N1)(1|(3 + 4)|2)(2N4) 
= (5|^ 2 l(2 + 3 + 4)2^ 2 |4) 

(5.3) 

= l\ (5|(3 + 4)|2)(2|^ 2 |4) + i\ (51) [\\i 2 {2 + 3)|4] + (51) [24] (1 + |£ 2 (2 + 3 + 4)£ 2 |2+) 

= t\ (5|(3 + 4)|2><2|* 2 |4> + £l (51) [1\£ 2 (2 + 3)|4] + (51) [24] t 234 (2\£ 2 \1) . 
The first two terms in eq. (5.3), after multiplication by (5|£ 2 |4), are triangle integrals with two 
external masses and two powers of the loop momenta in the numerator. They can be handled 
straightforwardly by Feynman parametrization. The third term is still a quadratic box integral 
(i.e. two loop momenta in the numerator), but it can be reduced further, using a second factor of 
(1|(3 + 4) 1 2) and a Fierz rearrangement (since £ 2 = 0), 

(1|(3 + 4)|2)(5|£ 2 |4)(2|£ 2 |1) = (5|£ 2 |1)(1|(3 + 4)|2)(2|£ 2 |4), (5.4) 
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and then following exactly the same steps as in eq. (5.3). 

Now we have five terms, only one of which is a box integral, and that one is merely linear in 
the loop momenta, 

I\ K 2 N1)] = (l|(3 + 4)|2) (terms lacking a t 234 cut). (5.5) 

Finally, we perform the four triangle integrals, dropping terms with no £234 branch cut, and assemble 
the pieces, thus obtaining all terms in eqs. (10.10) and (10.11) for Aq c (1^, , 3~, 4 + , 5J , 6+) that 
have logarithms or dilogarithms with £234 in the argument. (Terms representing possible scalar 
one- and two-mass triangle contributions are proportional to (— £234) -e /e 2 an d therefore can be 
inferred from the known structure of the poles in e [7,47]. This can be used as a check on the cut 
calculation, or to save labor.) 

The second cut we consider is the cut in the £412 channel of the amplitude Ag c (l+ , 2 + , 3^~, 4~ , 5^ , 6+). 
In this case leg 4 has to be adjacent to leg 1, as in the integral 1 5 in fig. 7. For this configuration we 
label the propagator momenta around the loop, starting just clockwise of the 5-6 lepton pair, by 
£' 5 , £1, £2, £3, ti- Now the required product of tree amplitudes, to be integrated over the two-body 
phase space for £' 5 + (—£3), is 

P 4 i2 = 4 ree ((4) g ",(-4) s ,3r,5 e r,6+)x4 ree ((-4)+ 4",1+ 2+,(4) s ) 

= [l 3 6] 2 (l 3 4) 2 (5-6) 

[43] [56] X (12) (2£ 3 ) (£ 3 £> 5 ) ' 

Fortunately, this cut can be obtained from the first one we evaluated. If we multiply eq. (5.2) by 
[34] / [24], and then perform the following 'flip' (pairwise permutation plus complex conjugation of 
all spinor products) 

£\ <-> £ 4 , £ 2 <-> £ 3 , £5 <-> £' 5 , 1 <-> 3, 5 <-> 6, {ab) <-> [ab] , (5-7) 
we recover eq. (5.6). 

This second example illustrates the principle of recycling cuts. One can save a lot of effort by 
identifying different cuts that are actually the same up to permutations and trivial overall factors. 
In many cases, it may not be possible to obtain an entire cut in this way, but portions of it may 
be recyclable. The 'master functions' defined in section 7, which enter three different amplitudes, 
are the most complicated contributions we have been able to recycle, but there are several other 
instances as well. 

5.2 Rational Function Reconstruction 

In Section 3.2 we described the general factorization properties of amplitudes, and how that 
information can be used to determine the rational function parts of amplitudes. As an example, 
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we explicitly construct the rational function terms proportional to n/ in the leading-color helicity 
amplitude, which are given by Ag 2+, 3 + , 4^~) in eq. (8.2), and outline the construction of the 
scalar (non-cut-constructible) piece (V sc , F sc ) for the same helicity configuration, eq. (8.7). 

The first step is to account for possible spurious singularities in the rational function terms. 
These singularities can always be identified after all cuts have been calculated, because they have 
to cancel against terms containing logarithms and dilogarithms. For A$ 2 + , 3 + , 4^), the only 
possible cut, that in the S23 channel, vanishes identically, because the tree- level q~q + g + g + and 
ssg + g + amplitudes vanish for massless quarks and scalars. So this term is purely a rational function, 
and therefore can have no spurious singularities, and we can immediately focus on the physical 
(multi-particle and collinear) singularities. 

It is usually convenient to match the multi-particle behavior — here, the limits £123,^234 — * 
- before attacking the collinear poles. The residue of a Ujk pole in a one-loop six-point amplitude 
receives three contributions:''" C\ from a one-loop four-point amplitude multiplying a tree-level four- 
point amplitude; C2 where the loop amplitude is replaced by the corresponding tree amplitude, and 
the tree by the loop; and C3 is associated with a loop correction to the intermediate propagator, 
multiplied by the two tree amplitudes. In the case of e + e~ — > four partons, two of the three 
contributions are easy to describe in general, because one of the two four-point amplitudes is the 
relatively simple e + e~ — > qq process. Let us define the C2 contribution to be where the e + e~ — > qq 
amplitude is a loop amplitude. Then for the n/-dependent term, C 2 and C3 both vanish. Also, for 
the scalar (sc) part of an amplitude, C2 = \ x lim^4g ree , and C3 = (—1) x lim^4g ree , where limylg ree 
is the appropriate tijk — ► limit of the six-point tree amplitude. 

For the present (1+, 2+, 3 + , 4^") helicity configuration, the tree amplitude (8.4) has no tijk 
poles, hence only the C\ contribution survives. Using a calculation of the n/ terms in the qqgg 
one-loop four-point amplitudes [15,19], we have 



ic r s 56 ti2 3 3 s 23 (12)(23) 



«56 ^234 3 S 23 (P2) (23) 



Similarly, the limits of the scalar pieces are found to be 



^ c (l+,2+,3+,4 g -) (45) [6P] 1 



x x 



i Cr «56 *123 

- [16] <5P) 1 

~ — X X 

S56 ^234 



11 £12 

2 3s 23 



(P|l|3) 



(12) (23) 

1 1 s 34 



(5.8) 



(P = 1 + 2 + 3), astiss^O, 



^'l P l' i; (p = 2 + 3 + 4), ast 23 4 ^0. 



2 3s 23 J <P2) (23) 

(5.9) 



t In the presence of infrared divergences, such as those due to virtual gluons, the one-loop factorization is a bit 
more subtle, but still has a universal form [30]. 
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Next we need to modify the terms in eq. (5.8) to improve their collinear limits, while preserving 
their tijk limits. For example, we can improve the k 2 || k 3 limit of the £123 limit in (5.8) as follows: 

1 (45) [6|(l + 2 + 3)l|3]si2 _ 1 <45) [6|(1 + 2 + 3)1|3] [12] 

(5.10) 



3 (12) (23) 523*56*123 



(23) 2 [23] s 56 t 



123 



l(45)[6|(l + 2 + 3)2|3][12] 



_ 1 (45) [6| (l + 2 + 3)2| 1] 



3 (23) 2 [23] 556*123 ' 3 (23)^56*123 

where nonsingular terms in the £123 — > limit are represented by '• • •'. A similar manipulation of 



the £234 limit, using also s 2 p = — S34 in that limit, gives 

l [16](5|(2 + 3 + 4)|2)(4|2|3) | _ 1 [16] (5|(2 + 3 + 4)2|4) | 

3 (23) S 23 556*234 3 (23) 2 S 56 t 2 34 

Thus a first guess for 2 + , 3 + , 4^)/(i c r ) would be 

1 1 r<45)[6|(l + 2 + 3)2|l] [16] (5|(2 + 3 + 4)2|4) 



(5.11) 



guess 



3 (23) 2 s 56 



!-123 



*234 



(5.12) 



It might appear that the k 2 \\ k 3 and k§ \\ k§ collinear limits of (5.12) both need to be improved, 
since the expected limits are ~ 1/ ^/SiJ, not 1 / Sij . But in fact there is a cancellation between the 
two terms in eq. (5.12) in both limits. In the k 2 \\ k 3 limit, using the identity 



S12 *234 - S24 *123 = S12 S34 - S24 «13 + «23 «14 + 0(s 23 ) = tr[2341] + C(s 23 ) 

= (23) [34] (41) [12] + [23] (34) [41] (12) + 0(s 23 ), 
and letting P = 2 + 3, we have 

"(45) [61] 



(5.13) 



F s -l—L 

^ guess Q /rlo \ 2 



3 (23) 2 S56 



*123*234 



((23) [34] (41) [12] + [23] (34) [41] (12)) 



_ (23) (45) [63] [21] + [16] (53) (24) 



1 v*a 


-z) 


3 (23) 2 






-z) 


3 (23) 2 


■S56 


1^(1 





~ (23) (45) 



*123 *234 

(4|1|6) + (4|P|6) 



(IP) (P4) 



+ [23] [16] 



(5|4|1) + (5|P|1) 
[IP] [P4] 



(23) + [23] [16]2(56) 



(23) 



(IP) (P4) 

~ (45) 2 
(IP) (P4) (56) 



+ 



[IP] [PA] 
1 y/z(l - z) [23] 



[16]' 



(23)' 



[IP] [P4] [56] 



k 2 || k 3 . 

(5.14) 

Not only does the leading 1 /|«23 1 term cancel in eq. (5.14), but the next terms in the expansion 
give precisely the desired limit, corresponding to the ny-dependent piece of the 



Splid- loop (2 + ,3+) x 4 ree (l+P ± ,4,-,5 e -,6+) 



(5.15) 
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terms in eq. (3.3). The loop splitting functions are given in ref. [22] and the five-point e + e~ — ► qqg 
amplitudes in appendix IV. Notice that there are no nf terms in Al~ loop , hence no Split tree xAl~ loop 
terms contribute here. Similar manipulations show that the k$ \\ k& limit of P| uess is also precisely 
correct, 

F> 1-^,1 [23] (24) (34) , z 1 [12] [13] 

'guess ~ [56] 3 (23)2 (4p) (pi) + (56) 3 [4p] [pi] (23) , *5 || *«. (5.16) 

Here there are only Split tree xAg _loop terms, with the relevant loop amplitudes obtainable from 
ref. [19]. All singular collinear limits of P| uess have now been verified, so we expect the result to 
be correct as is. A simple identity shows that indeed P| uess = Aq(1+ , 2+, 3 + , 4^")/(i cr), the result 
given in eq. (8.2). 

A slightly more complicated example, which we only summarize here, is the rational function 
terms in the scalar piece (V sc , F sc ) for the same helicity configuration, ^4g c (l^", 2 + , 3 + , 4^"), as 
given in eq. (8.7). By comparing the limits (5.8) and (5.9), and similarly the collinear limits, 
it is easy to see that there must be a term of the form Aq in Aq c . As for the remaining terms, 
the first step is again to account for possible spurious singularities. Assuming that all the cuts 
have previously been calculated, the ln( )/( s ij ~ ^234) 2 terms contained in the Li functions in 
eq. (8.7) are already known. But the Li functions are designed to cancel the spurious (sjj — £234) 
behavior between logarithms and rational functions. Hence by completing the ln( _ f ^ ) / (sij — 1234) 2 
terms into Li( _^' : ' 4 )/^234 functions, we ensure that the remaining rational function terms are free 
of l/(sij — £234) poles. Since there are no other non-rational-function terms in V sc or F sc , there 
are no other sources of spurious poles, and we next turn to the physical singularities. 

Using the t-channel limits (5.9), and manipulations similar to those leading to eq. (5.12), we 
arrive at an ansatz for what has to be added to the L 1 terms, 



F sc = - 

guess 2 



(45) [6|(2 + 3)1|3] + [23] [16] (5|(2 + 4)|3) (5|2|3) 2 , <4|(I 



(12) (23) t 123 s 56 (23) [34] t 234 s 56 (12) (23) [34] (56) t 234 (12) (23) (34) [56] t 234 

(5.17) " 

The first two terms reproduce (5.9), while the last two terms are needed to cancel off £234 poles in 
the Li terms. 

This time, however, the k§ \\ limit, 



F SC r I 

x guess 2 



(45) [64] (41) [13] [23] [16] (51) [13] 
(12) (23) i 123 s 5 6 (23) [34] t 234 s 56 



Vgj ~ g [13] 
2 (23) s 56 



<4P) [P4] (41) + [23] [IP] (PI) 



(12) s 4P [34] sip 



(5.18) 



Vgl^g [3P](P1)[13] 
2 (12) (23) [34] s 56 ' 
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still has to be improved, for example by adding 

SFr = l [36] (5|(2 + 4)|3) 

6tl 2 (12) (23) [34] s 56 l "- iyj 

to -Fgu ess - At this stage, the polynomial ansatz formed by the Li terms, F^ ess and 5F( C has only 
collinear singularities. One can now systematically add additional terms to match the 
known ki \\ k 2 1 1 ^3, ks \\ & 4 , and k 5 \\ k§ limits. (One can also rewrite the answer in a form 
where the I/^/s^q behavior is manifest, as in eq. (8.7).) 

Other helicity amplitudes may possess more spurious singularities and/or physical singularities. 
For these cases the determination of the rational function terms becomes somewhat more involved, 
but the basic principles remain as illustrated above. 

5.3 Numerical Simplification 

In order to illustrate the numerical simplification technique outlined in section 4.3, we consider 
the particular example of the coefficient c 3m of the three-mass triangle integral /| m (si2, S34, s^q) in 
the axial- vector fermion-loop contribution yl ax (l+, 2^~ , 3 + , 4~), eq. (11.9). Note that the three-mass 
triangle contribution is contained in the common function C ax defined in eq. (11.8), plus its image 
under the symmetry operation flip 2 defined in eq. (6.8). Thus we determine simultaneously the 
corresponding contribution to A ax (l+, 2^~, 3~ , 4 + ), eq. (11.10). 

The first step in the technique is to write down terms reproducing all of the spurious and 
physical singularities in the various channels. One subtlety is that a term reproducing a singularity 
in one channel may be too singular, or otherwise have the wrong type of singularity, in another 
channel. In this case the term will have to be 'improved'. 

In the present example of c 3m , we first write down a term cf m which reproduces the known 
(3|(1 + 2)|4) — > behavior of this coefficient. To do this we use eq. (4.22) and the coefficient of the 
hard- two- mass box function in C ax (which we assume has already been obtained via the £123 cut), 
to get 



l (2|(l + 3)|4) 2 (3|(l + 2)|6) 2 -(23)Wtr2 3 
tess 2 (12)[56](3|(1 + 2)|4) 4 

/W1 M 1\\A\/A\(1 J- 9^U ^ 3 _ tl24 ) Sl2S56 



Ci ' & "™ 2 (12) [5t)J (31(1 + 2j|4)' ± 

(3|(1 + 2)|4)(4|(1 + 

t 2 A 



123^3 

(2|(1 + 3)|6) [12] (23) [46] (56) (t 123 - t 124 ) (4|(1 + 2)|3) 
ii 23 (3|(1 + 2)|4) 2 A 3 

In the second line of eq. (5.20) we used the spinor identity 



+ 



(2|(1 + 3)|4)(3|(1 + 2)|6) = (3|(1 + 2)|4)(2|(1 + 3)|6) - (23) [46] t 123 , (5.21) 
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and dropped all but the leading terms as (3|(1 + 2)|4) — > 0. As it stands, c 3 ™ uess contains a pole in 
ti23, but an identity similar to (5.21) removes the pole while preserving the leading (3|(1 + 2)|4) — > 
behavior, and so we take the first singular term to be 

am _ [14] (35) (t 123 -t 124 )(4|(l + 2)|3)(2|(l + 3)|6) 
1 " (3|(1 + 2)|4)2 A 3 • 

Actually c 3m should contain a pole in £123 (but with a different structure than that found in 
cq. (5.20)); it has to cancel the pole in the explicit formula for the hard two-mass box function, 
eq. (II. 3). Thus a term reproducing that pole is given by 

c3m ,1 *12*56 (2|(1+3)|6) 2 

2 ' guess 2 t 123 (12)[56](3|(1 + 2)|4) 2 ^ ' { } 

where we used identity (5.21) again and dropped nonsingular terms as £123 — > 0. We still have to 
remove the leading singularity as (3|(1 + 2)|4) — > in this term, which we can do using 

(1+|2(1 + 3)6|5 + ) = -(1+|3(1 + 2 + 3)6|5+) + • • • = (1 + |3(1 + 2)4|5+) + • • • , (5.24) 

as £123 — > 0, thus obtaining 

3m _ 1 [13] (45) (2|(1 + 3)|6) 

C2 " 2 t 123 (3|(l + 2)|4) (5 ' 25) 
as our second inferred singular term in c 3m . 

As mentioned above, the leading singular terms as A 3 — > require an explicit calculation of 
the Si2-S34-S56 triple cut, but fortunately only the leading loop-momentum terms in this cut have 
to be retained, yielding 

3m £34 «5|2|1)* 12 - (5|6|1)<5 56 )(4|(1 + 2)|3)(2|(1 +3)|6) 
° 3 "" 3 (3|(1 + 2)|4) A 2 • (5 - 26) 

The remaining terms can have at most a I/A3 singularity. Furthermore, there are no collinear 
singularities (ki \\ kj) in the three-mass triangle coefficient. (This feature is general, as long as the 

2 l 2mh 

adjacent two-mass box function Ls_ 1 is used, and not Ls_ 1 .) 

At this stage we have identified enough of the singularities of c 3m to write an ansatz for the 
remaining, less singular terms. The unknown coefficients in the ansatz can then be determined 
by comparing it to a numerical evaluation of c 3m at a number of phase-space points, equal to the 
number of unknown coefficients. The full three-mass triangle coefficient is symmetric under both 
the operations flip 2 and flip 3 , defined in eqs. (6.8) and (6.9), and the full singular terms include 
also the images under flip 2 of the above terms. (They cancel against terms containing the other 
box function, Ls 2 _™ /l (s 3 4, £124; S12, s 5e ).) Thus we may write 

3m 3m i 3m . 3m i Pi , ( 4 I(1 + 2 )I 3 ) , fl; / c -ov^ 

c _ Cl +c 2 +c 3 + — +P2 _____+ ttlp2 , (5 .27) 
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where flip 2 is to be applied to all preceding terms. In eq. (5.27) we have explicitly identified 
all denominator factors; i.e., p\ and p 2 can only be linear combinations of products of spinor 
strings. They are further restricted by the observation that any helicity amplitude can be assigned 
a definite phase weight — namely, an integer n, for each external leg i, which is equal to twice 
the helicity of that leg. For ^4 ax (l+, , 3 + , 4~, 5J , 6+), and hence for C ax , the phase weight is 
{1, —1, 2, —2, —1, 1}, the six entries corresponding to each of the six legs. The phase weight of an 
expression A may also be calculated by making the substitution (ij) — ► A~ 1 A~ 1 (ij), [ij] — ► AjAj [ij], 
and reading the rii off of the resulting transformation A — > (FJ i A"* ) x A. 

Since A 3 carries no phase weight, the phase weight of p\ must match that of C ax , namely 
{1,-1,2,— 2,— 1,1}. In terms of spinor string terminations, it must have the form 

[1, (2, ([3) 2 , «4) 2 , (5, [6 . (5.28) 

On the other hand, dimensional analysis implies that the mass dimension of p 1 is 4. (The mass 
dimension of C ax , like that of any six-point amplitude, is —2; but /| m also has dimension —2, 
and A 3 has dimension 4.) Thus there are exactly four spinor products ((ij) , [kl]) in p 1 , all of 
whose arguments are accounted for in eq. (5.28). Since (ii) = [ii] = 0, we see that p\ has just one 
independent term, 

pi = a [13] [36] (24) (45) , (5.29) 

where a is a constant. 

Similarly, p 2 also has dimension 4, but its phase weight is {1, —1, 0, 0, —1, 1}, corresponding to 
the spinor string terminations 

[1, (2, (5, [6 . (5.30) 

Equation (5.30) allows terms of the form (5|i|6)(2|j|l), (5|i|l)(2|j|6), (52) [l|ij|6], (5|ij|2) [16] and 
(52) [16] Sij, where i and j are arbitrary legs. Not all such terms are independent. A Fierz identity 
lets us eliminate all terms of the first type in favor of terms of the second and fourth types. The 
symmetry flip 2 relates terms of the third and fourth types, while combinations of these terms 
that are symmetrized in i <-> j are already included in the fifth type of terms. Using momentum 
conservation and the two flip symmetries, we can reduce the number of independent terms in p 2 to 
just seven, 

p 2 = ai (5|3|1)(2|3|6) + a 2 (5|3|1)(2|4|6) + a 3 (5|31|2) [16] 

(5.31) 

+ (52) [16] (a 4 S12 + a 5 s 23 + a 6 s 24 + a 7 s 25 ) + fhp 3 • 
Thus it suffices to evaluate c 3m numerically at eight phase-space points, and solve the resulting 
linear equations for ao, cq, . . . , 07, which should be simple rational numbers. We find 
p 1 = - [13] (45) (24) [36] , 

P2 = -|(<5|2|1> <2|1|6> + (5|6|1)(2|5|6) - (5|3|1)(2|4|6) - (5|4|1)(2|3|6)) . ^ 
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(Evaluation at additional phase-space points may be used to confirm this answer.) If we combine 
these terms with the singular terms cf m according to eq. (5.27), we obtain the coefficient of the 
three-mass triangle given in eq. (11.9), via the function C ax defined in eq. (11.8). 

This numerical simplification procedure may be somewhat more involved for other coefficient 
functions, if there are more spurious and physical singularities to account for, but the principle is 
the same. In practice we found it convenient to obtain an analytic, though often very complicated, 
representation of the desired answer, before attempting to simplify it numerically; of course, an 
analytic representation is not strictly necessary. It is also not necessary to eliminate all the linear 
dependences between terms in the ansatz before trying to solve for the unknowns. 



6. General Form of Primitive Amplitudes 



The simple structure of the poles in e of the primitive amplitudes [7,47] permits us to decompose 
them further into divergent (V) and finite (F) pieces, 



Al- loop = c r 



A% ee V + i F 



(6.1) 



where the A\ loop are any of the primitive amplitudes and the prefactor is 



cr 



1 r(l + e)r 2 (l-e) 



(47T) 



2-e 



r(l - 2e) 



(6.2) 



The tree amplitudes [5] are denoted by A§ ee . (The Aq^'* amplitudes are rather simple so we will 
not bother with this additional decomposition for them.) 

The amplitudes we present are bare ones, i.e., no ultraviolet subtraction has been performed. 
To obtain the renormalized amplitudes in an MS-type subtraction scheme, one should subtract the 
quantity 



c r N c g 2 



2 n f 
3~K 



3N r 



A, 



tree 



(6.3) 



1 /ll 

from the amplitudes (2.9) and (2.11). 

We quote the results in the four-dimensional helicity (FDH) scheme [43,44], since this scheme 
is convenient for performing computations in the helicity formulation. The conversion between the 
various schemes is discussed in refs. [43,15]. (The more conventional regularization schemes alter 
the number of gluon polarizations and are therefore not natural when using a helicity basis.) The 
e + e~ — > qqgg amplitudes may be converted to the 't Hooft-Veltman (HV) scheme by adding the 
quantity 



\c T N c g 2 



1 



1 - — I ^e" 6 



(6.4) 
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to the amplitude (2.9) and changing the coupling constant from the non-standard ajj^ to the 
standard ctj^g. Similarly, the e + e~ — > qqQQ amplitudes may be converted by adding the quantity 

-^ c <7 2 0-^)^ 6 ree (6-5) 

to the amplitude (2.11) and making the same coupling constant conversion. The conversion of 
the HV scheme to the conventional dimensional regularization (CDR) scheme is accomplished by 
accounting for the fact that in the HV scheme observed gluons (at the partonic level) are in four- 
dimensions but in the CDR scheme they are in (4 — 2e) dimensions. This conversion is rather simple 
as it involves only the coefficients of the poles in e which are proportional to the tree amplitudes. 
In the final matrix elements squared (e.g., eq. (2.12)) one simply replaces all terms originating from 
the singular parts of eqs. (6.1) and (6.3) with their values in the CDR scheme. (Further details may 
be found in ref. [15].) Using this conversion recipe, our results (integrated over lepton orientation) 
agree numerically [48] with those reported in refs. [13,14]. 

In order to present the amplitudes compactly, we define, in addition to the previously-defined 
spinor-strings, the following combinations of kinematic variables, related to the three- mass triangle 
integrals that appear, 

<5l2 = Sl2 — S34 — S56 , 834 = S34 — S56 — Sl2 , #56 = S56 — S12 — S34 , 

^3 = s\ 2 + £34 + Sg 6 - 2si 2 s 3 4 - 2s 34 s 56 - 2S 56 S12 , 

(6.6) 

^3 = S?4 + s 23 + s 56 ~ 2si 4 S 23 ~ 2s 23 S56 ~ 2s 56 Si 4 . 

Certain 'flip' symmetries relate either various terms, or various amplitudes. For later conve- 
nience, we collect the definitions of these symmetry operations here: 

(6.7) 

(6.8) 
(6.9) 
(6.10) 

In general, {ab} <-> [ab] denotes complex conjugation of all spinor products, including the various 
strings defined in eq. (2.2), <-> (l\j\i), (i\(l + m)\j) <-> [i\(l + m)\j], and so forth. We also define 

the exchange operations 

exchi 6 25 : 1 6 , 2^5, 

(6.11) 

exch 34 : 3^4. 



fliPi : 


1 ^4, 


2 


<-> 


3, 


5 <-> 6 , 


{ab} <-> 


[ab] 


flip 2 : 


1^2, 


3 


<— > 


4, 


5^6, 


(ab) <-> 


[ab] 


flip 3 : 


1 ^5, 


2 


<-> 


6, 


3^4, 


(ab) <-> 


[ab] 


flip 4 


: 1 3, 




5 ■ 


i ^ 6 , 


{ab} * 


-> [ab] , 
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7. Master Functions 



As mentioned in section 5.1, many of the cuts (or portions of cuts) for different amplitudes are 
related to each other by simple permutations and overall spinor-product prefactors. In particular, 
one 'master function' M\{1, 2, 3, 4) enters the the cut-constructible part of one of the leading-color 
primitive amplitudes, as well as two of the subleading-color ones. Two additional 'master functions', 
M2(l, 2, 3, 4) and M^fl, 2, 3, 4), enter the scalar parts of the three subleading-color amplitudes with 
opposite gluon helicities. We have decomposed the latter two master functions further, extracting 
M 2a and M 3a , respectively, because these pieces appear separately as well. 



The cut-constructible master function is given by 



Mr (1,2, 3, 4) = 

(1|(3 + 4)|2)(3|(1 + 2)|4)A 3 ( 2 {U) < 5 I 2 ' 6 ^ 123 Sl2 + 556 ^ + 2 < 12 > < 15 > < 4 l 5 l 6 >(t 14 ] ^ " 2 I 1 ! 23 ! 4 ]) 
+ s 56 (24) ((15) [46] <5 56 -2(l|2|6)(5|3|4)))/3 3m ( Sl2 , S3 4,5 56 ) 

[13] 



+ 



2 



■ «2|1(3 + 4)| 5) - (2|(3 + 4)(1 + 2)|5>) ((15) i 124 + <1|26|6» 



(1|(3 + 4)|2)(3|(1 + 2)|4)A 3 L(56) 

[J ° ] 2 [) ((*i4-S2 3 )«l|2|6>tf 12 - (1|5|6)<5 56 ) + (1|(3 + 4)|2)((2|1|6)<5 12 - (2|5|6)<5, 



(34) [56] 

+ 2(2|l|3)(5|(l+2)|6)( Sl4 -s 2 3) 



.-.(> 



In 



(7.1) 
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The scalar master functions are given by 



M 2o (l,2,3,4) 



1 [12] J 6 (12) (3|(1 + 2)|4)(5|(1 + 2)|6)([13] <5 56 - 2 [1|24|3]) 



2 [23] (4|(1 + 2)|3)\ A§ 



1 

+ A 3 



(2|1|4)(35) (5|(1 + 2)|3) (12) (45) / [13] (45) (3| (1 + 2)\A)5 12 M4U hiooMi^ 
(56) (34R56)l <4|(1 + 2)|3) (35) (2 - [1|23|4]) j 

^|^((5|(1 + 2)|4) - (5|6|4)) + ||(2 ((2|3|6)«5 34 - (2|5|6)<5 56 ) + (3|4|6)(2|(1 + 4)|3)) 

^w+w^ " 4 (12) [36] ( (35) [14] + + ^ « 2 ' 46 ' 5 > - < 2 ' 13 ' 5 » 



(5|(1 + 2)|3) ((25) (34) -(23) (45)) ) / x 
(34)(56)(4|(1 + 2)|3) J v"W 



1 1 / _ 6 [12](5|(1 + 2)|6)[4|3(1 + 2)|4]((24)«5 56 - 2 (2|13|4)) 



2 [23] <4|(1 + 2)|3> A 3 I A 3 

[13] , x/ (5|34|5) [6|34|6] 



(4|(l+2)|3) v MV " d 1Z " ; a/ V (56) [56] 

+ lif ( 2 ((5|4|1) 534 " (5|6|1) * m) + (4|(2 + 3)|1)(5|(1 + 2)|4) ) 

- {H}( 2 <3|2|1)([36] *w " 2 PI 45 ! 6 !) + m + 3)43(1 + 2)|6] - 3 [6|43|1] <5 34 ) } In 



(7.2) 



M 2 (l,2,3,4) 



M H 2 341 l (4|(2 + 3)|1)(5|(1 + 2)|3) 2 2mfe 
M2a(1 ' 2 ' 3 ' 4) + [23](56)(4|(1 + 2)|3) 3 ^ (S34 ' * 123 ' Sl2 ' ^ 

[12] f (12) [34] (3|(1 + 2)|4)(5|(1 + 2)|6) , , 
+ [23] (4|(1 + 2)|3) \ 3 A| ( (4|2|1) 512 " (4|3|1) ^ 



1 

A 



(2| L3( ! + 2)4|5) [36] ^ _ ^ + (2g) [4fi] ^ _ (g|6|g) ^ 



(4|(1 + 2)|3) 

+ (3|(1 + 2)|4) ((2|4|3)(5|(1 + 2)|6) + 3 (2|1|3)(5|4|6)) 
- (2|1|4) (35) ([6|4(1 + 2)|3] + 2 [6|54|3] + [36] <5 12 ) 



(2|3(1 + 2)|5)[46H f3m , , , x 

H 7 M 3 (Sl2, «34, «56j 

C123 J 



[46] <5|4|1> L H^SJ (3|2|1)(5|(1+2)|3) 2 L o (5|4|1)(5|(1 + 2)|3)t 123 L o 



[23] <4|(1 + 2)|3> t 123 [23](56)(4|(1 + 2)|3) 2 s 12 [23] (56) (4|(1 + 2)|3) 2 s 56 

1 (4|2|1)<5 12 -(4|3|1)<5 34 /(34)[46] 2 (35) 2 [34K (35) [46] ([13] <5 56 - 2 [1|24|3]) 



+ 



2 [23](34)(4|(1 + 2)|3)A 3 V [56] (56) 7 [23] (4|(1 + 2)|3) A 3 

(7.3) 
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M 3a (l,2,3,4) 
1 



[24] 



2 [23]S34<1|(3 + 4)|2> 



6 £ ^ i (ti34-t 23 4)(3|(l + 2)|4)(5|(l + 2)|6) 



+ (m + 2 )l 6 ) t 14 ] I 26 ] + <3|(2 + 4)|6) [12] [46] - (13) [24] ([16] 2 + [26] 2 



+ 



- 2 



[56] A 3 

[12] 6 12 
(56) A 3 

(12) (35 



(5|2|4)((12) (35) - (23) (15)) - (15) (25) (3|(1 + 2)|4) + { -^^^^^ 



A, 



([16] [24] + [14] [26])<5 56 



(34) [46] 2 , (35) 2 [34]\ s 12 t 234 



<3|2|6> [46] (13) [24] / (5|12|5) 
[56] " (l|(3 + 4)|2)l (56) 



[56] 
[6|12|6] 
[56] 



+ 



In 



(56) 



A, 



-S56 



1 

+ ^ 



[24] 



5 5 6 (tl34-t234)<3|(l+2)|4)<5|(l+2)|6) 



2 [23] (1|(3 + 4)|2)A 3 L 

/(3|2|6)[46] (5|1|4)(35) 



+ 2. 



[56] (56) ' -<56><1|(3 + 4)|2> 

(13)[24] S56 \ /(5|12|5) [6J12J6K 

[56] J 



+ ((3|(l + 2)|4) + 4 



(1|(3 + 4)|2)M (56) 



™_((3|(2 + 4)1|5)-(3|26|5))). 34 
s 12 (15) [24] (3|5|6) 



+ 8- 



(1|(3 + 4)|2) 



2 (3|(1 - 2)|4) (5|(1 + 2)|6) + 4 ([4|12|6] (35) + (5|21|3) [46]) 



In 



-S34 

-S56 



(7.4) 



M 3 (l,2,3,4) 



M 3a (l, 2, 3, 4) (15) [24] t 234 2 h f , 

K ' [23] [34] (56) (1 1 (3 + 4) 1 2) 3 _1 V ' 

[24] f 3 Sl2 <5 12 (* 134 - * 2 3 4 )<3|(1 + 2)|4><5|(1 + 2)|6> + 1 <3|(1 + 2)|4)<5|(1 + 2)|6> 



[23](1|(3 + 4)|2)A 3 12 A 3 2 

" (13)(1 (l ) |[flT 2 ) +4)|1) ([2|1(3 + 4) ' 6] " [2 ' (3 + 4)(1 + 2) ' 6]) 
+ (2|1|4)((3|(1 + 4)|2)(5|(1 + 2)|6) + (5|1|2)(3|4|6)) 



+ (15) (23) ([24] ([1|2(3 + 4)|6] - [1|(3 + 4)(1 + 2)|6]) - [12] [4|(1 + 2)(3 + 4)|6] J }/ 3 3m (s 12 , s M , s 56 ) 



1 (15)' [14] [24] / [24] /-. 56 x [u1 Ll (^S) 

^2 [23][34](56)(1|(3 + 4)|2) V <1|(3 + 4)|2) °V-t 234 ^ L J 



^234 



_ 1 (23) [24] (5|(3 + 4)|2) // (5|(2 + 3)|4) _ (15) [24] \ L o(^f) L i(^gQ \ 
2 [23](56)(1|(3 + 4)|2)VV *234 <1 1 (3 -h 4) 1 2> ^ , 34 ) 

;! ' ! 2;>) «3|1|2) 6 12 - (3|4|2)<5 34 ) - [16] ((13) ,5 56 - 2 <1|24|3») 



, 1 [24] 
+ 2 [23] (1|(3 + 4)|2)A 3 



(34)V(56) 



+ [46] f[16] 



[34]V[56] 



«1|2|4) 5 12 - (1|3|4) <5 34 ) + (25) ([24] 5 56 - 2 [2|13|4])) 



(7.5) 
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8. Results for Primitive Amplitudes: A 6 (l q , 2, 3, 4$) and A s 6 ,f,t (l q , 2, 3, 4,) 

In this section we present the independent gggf/ primitive amplitudes where in the parent 
diagrams neither of external gluons are attached to the external fermion lines, i.e. the helicity 
configurations 

1+ 2+,3+,4r,5=-,6+; 1+ 2+, 3", 4r, 5=", 6+; 1+ 2", 3+, 4r, 5=", 6+ (8.1) 

We suppress the lepton labels (5 J, 6+) below. Representative parent diagrams are given in figs, la 
and b. 

8.1 The Primitive Amplitudes: A 8 /' 1 (1„ 2, 3, 4 ? ) 

By far the simplest of the primitive amplitudes are those proportional to the number of scalars 
n s or fermions rif. This simplicity follows from the fact that only two Feynman diagrams contribute, 
each with a triangle integral. (See fig. la.) It may also be understood in terms of the rather simple 
cut and factorization properties that must be satisfied. 

The results for these primitive amplitude are 

,4|(1+ , 2+ , 3" , 4r ) = 4(1+ , 2+ 3" 4r) = , 
A s 6 (l+, 2-,3+,4r) = 4(l+ 2" 3+ 4r)=0, 

4(1+ 2+,3+,4r)= , ( 8 - 2 ) 



A*(l+,2+,3+,4 g - 



3 (23} 2 s 56 



(45) [6|(1 + 2)3|1] | [16] <5|(4 + 2)3|4> 

*123 *234 



Although it is not manifest, A|(l + , 2 + , 3 + , 4^ ) is antisymmetric in the exchange of 2 and 3, as 
required by charge conjugation invariance. This fact accounts for the absence of n/ or n s terms in 
A 6;3 in eq. (2.13). 

The contribution of a virtual top quark, through bubble and triangle graphs, is simply related 



to the above function Ag, 



A*(l+,2+,3-,4r)=A^(l+,2-,3+,4r) = 0, 

4(l g + ,2+,3+,4r) = l£|^(l+,2+,3+,4r), (8 ' 3) 

neglecting l/m| corrections. 

8.2 The Helicity Configuration q + g + g + q~ 

The simplest of the non-n Sj / helicity amplitudes is the one where both gluons have the same 
helicity, A 6 (l + , 2 + , 3 + , 4^). The tree amplitude in this case is 



4 tree _ ■ (fj) /on 

6 (12) (23) (34) (56) ' 1 ' ' 
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The contributions to the amplitude in terms of the decomposition (6.1) are 



V cc = 



-«12 



+ 



-S23 



+ 



-S34 



A 4 



-«56 



-4, 



(8.5) 



/Itree 
^6 



+ 2 



La-J^fu-,^*-) -Ls_if 

\ — *123 ' —*123 / \ 

(45)(5|2|3) L ° 



-S23 ~ $34 
"*234 ' — *234 



-S34 



(12) (23) (56) s 34 



+ 2 



^ — Ls?.™ e (ti23, *234; ^23, 5 56 ) 

(45) (5|1(2 + 3)|4) L o 



(12) (23) (34) (56) t 234 



V sc = — 



1 

2e 



-«56 



1 

2 ' 



_F SC = 



^tree 



1 

2 
+ 



fc 234 



1 / (4|(2 + 3)1|5) 

2 (45) 



+ 



_1 / <4|32|5) \ 2 L i 

(5|2|3) (54) 
(12) (23) i 234 (56) 

(4|(2 + 3)|6) 2 
(12) (23) (34) t 234 [56] ' (12) (23) t 123 t 234 [56] 



2 Li 



t 2 

t 234 



[23] (45) (5|(2 + 4)|3) 
(12) W234 (56) 

(4|2|3)[6|1(2+3)|6] 



+ 



+ — A s 6 (l+, 2+,3+,4r 
1 cr 



(8.6) 



(8.7) 



8.3 The Helicity Configuration q + g + g q 

The next simplest helicity configuration is Aq(1^ , 2 + , 3~, 4^~). Notice that this amplitude is 
symmetric under the 'flip]/ symmetry in eq. (6.7). The tree amplitude for this helicity configuration 



is 



^ ree = i 
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(12) [34] S23 S56 

An alternate form for A^, which has manifest behavior in the collinear limits k 2 \\ k 3 and k 5 || fc 6 , 
at the expense of more obscure behavior on the three-particle poles, is 



Al ree = i 
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The results for the cut-constructible and scalar pieces are 
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2 (12) [23] [34] [56] 

where 'flip 1 ' is to be applied to all preceding terms in the given expression. 

8.4 The Helicity Configuration q + g~ g + q~ 

The most complicated leading-color helicity configuration is Aq(1+ 2~ , 3 + , 4^~). In general, 
those amplitudes where the negative and positive helicities alternate around the loop are the most 
complicated ones. The symmetry flip 1 (6.7) holds for this helicity configuration as well. The tree 



amplitude in this case is 
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An alternate form with more manifest behavior in two-particle channels, but less manifest behavior 
in three-particle channels, is 
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The results for the cut-constructible pieces are 
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The results where the gluon in the loop is replaced by a scalar are 
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9. Results for Primitive Amplitudes: 45(1^,2,3^,4) 

In this section we present the independent A&{\ q , 2, 3g, 4) primitive amplitudes, corresponding 
to the helicity configurations 

l+,2+,3 g -,4+,5 e -,6+; 1+ 2+, 3J, 4" , 5 e ", 6+; (9.1) 
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note that the configuration 1+, 2~ , 3^~, 4+, 5^ , 6+ is obtained from 1+, 2 + , 3^~ , 4~ , 5jT, 6+ by the 
operation 'flip 4 ' defined in eq. (6.10). For these amplitudes one of the external gluons is cyclicly 
adjacent to the vector boson, as depicted in fig. lc. These primitive amplitudes only contribute to 
subleading-in-color terms. 

9.1 The Helicity Configuration q + g + q~ g + 

We now give the primitive amplitudes for Aq(1^ , 2 + , 3^~, 4 + ). These amplitudes are relatively 
simple because the three-mass triangle does not appear. The tree amplitude is 
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(9.4) 

The terms where a scalar replaces the gluon are 
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9.2 The Helicity Configuration q + g + q g 



We now present the primitive amplitude A 6 (l+, 2 + , 3^ , 4 ). The tree amplitude is 
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An alternate form is 
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The first form has manifest t-pole behavior while the second form has manifest behavior for fc 5 || k§. 
The cut-constructible contributions are 
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where the three-mass-triangle coefficients are 
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The contributions where a scalar replaces the gluon in the loop are 
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10. Results for Primitive Amplitudes: 2q, 3, 4) 

In this section we present the A e (l q , 2q, 3, 4) primitive amplitudes, corresponding to the helicity 
configurations 

l<j~> 2q , 3 + , 4 + , 5 § , 6^; l^~,2^,3 + ,4 ,5 E ,6^; 1^,2^,3 , 4 + , 5 g , 6|j~ ; 

again we shall suppress the lepton labels below. In the parent diagrams for these amplitudes, both 
external gluons are attached to the same fermion line as the vector boson, as depicted in fig. Id. 
These primitive amplitudes only contribute to subleading-in-color terms. 

10.1 The Helicity Configuration q + q~ g + g + 

Here we present the primitive amplitude yl 6 (l+, 2^~, 3 + , 4 + ). The tree amplitude is 
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(34) (41) (56) - 23 (23) (34) (41) (56) 

The contributions where a scalar replaces the gluon are 

1 



«56 



V s 



2e\-s 56 J 



+ 



2 ' 



(10.4) 



_p sc = — . 



(24) 2 (51) 2 



(23) (34) (4lT(56) 



Ls 2 ™ 6 (t\23, *234; -23 , -56) 



(23) (51)' 



-Ls_ 



-«23 



(23) (54)' 



+ 



(34) (41) (13)" (56) -V-*i*»'-ti M y (34) J (41) (56) 
((41)(52)-(24)(51))(5|4|3) Lo(^f) 1 (5|4|3) 2 (32) L 



Ls 2 ™ 6 (ti24, ^123! s 12, -56) 



(52)' 



(34) (41) 2 (56) 
2 

I ~*234 ) 

V -S56 ; 

) 



-23 



2 (34) (41) (56) 



t 2 

r 234 



2 (23) (34) (41) (56) 
1 (56)(2|1|6) 2 L i 



In 



+ 



+ 



(23) (34) {Aiy (56) 
(23) (5|(1 + 2)3|5) L °' 



(23) (34) (41) 



(24)(51) z (2|(3 + 4)|l) (52)(2|1|6) \ L o 



-^234 
-«56 



S56 



(12) 2 (5|(2 + 3)1|5) L o 



+ 



2 (23) (34) (41) s§ 6 (34) z (13) (56) s 12 (23) (41) z (13) (56) «23 

f (24) (5|4|6) _ (2 (41) (52) + (51) (24))(2|4|6) (52) (54) (2|(1 + 3)|4) \ L o(^f) 
V (34) 2 (41) (23) (34) (41) 2 (23) (34) (41) (56) J - 56 

(52) 2 ln (^W)_ (56)(2|4|6) 2 L i(^f) (24) (54) (5| (1 + 2)|4) L o (^) 



+ 



(23) (34) (41) (56) 

(24) (5|3|6) L ° 
(34) 2 (41) ~~ 

(51) (52) 



(23) (34) (41) 



<34>" (41) (56) 



-12 



-56 
/ 1 



(41) (13) (56) V (34) V-*"/ 



(23) [36] 2 (56) 
(34) (41) 
(12) 



s 2 
*56 



(23) (41) 



In 



+ 



1 



[34] (4|(1 + 3)|6)(2|(3 + 4)|6) | [34] (52) (5|1|3) | 



(52) (5|4|3) 



(34) (41) i 234 ti34 [56] 



(41) t 23 4ti34 (56) (34) (41) t 234 (56) 



+ 



(25)' 



(23) (34) (41) (56). 

(10.5) 



10.2 The Helicity Configuration q + q g g + 

Here we present the primitive amplitude A 6 (l+ , 2^ , 3~ , 4 + ). This amplitude satisfies the sym- 
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metry 'flip 2 ' defined in eq. (6.8). The tree amplitude is 



A tree 
A 6 



[41](13)(52)(3|(1 + 4)|6) , (23) [24] [61] (5 1 (2 + 3) 1 4) (5|(2 + 3)|4)(3|(1 + 4)|6) 



(41) S 3 4tl34S56 



+ 



[23] S34*234S56 



[23] (41) S34S56 



(10.6) 



An alternate form with manifest k§ \\ k§ behavior is 



a tree 
^6 



— I 



[41][24](52)(5|(2 + 3)|4> , (23) (13) [61] (3| (1 + 4)|6) , (3|(1 + 4)|6)(5|(2 + 3)|4) 



[23] [34] t 134 t234 (56) 



+ 



(34) (41) t 134 t 234 [56] 



+ 



[23] (41) i 134 i 2 34 



(10.7) 



The cut-constructible pieces are, 



V cc = 



-Sl2 



2 

e 



-S56 



(10.8) 
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(5|(2 + 3)|1) 2 /~2mh 

IjS 



[23] (56)(4|(1 + 2)|3)(4|(2 + 3)|1) 

+ 



1 (S14, *123l «23, S 5e ) + Ls_i (^^j 



[24] 3 (15) 2 t 234 [24](1|(2 + 3)|4) 2 (5|(3 + 4)|2) 2 

[23] [34] (56) (1|(3 + 4)|2) 3 [23] [34] (56) t 234 <1| (3 + 4)|2) 3 



(23) 2 [16] 2 (3|(2 + 4)|1) 



(34) [56] t 234 (2|(3 + 4)|1)(4|(2 + 3)|1) 



\ Ls 2 ^ h (s 12 ,t 2 34;S34,S 56 ) 



, [13] 2 (45) 2 ^ 23 (4|(2 + 3)|1)(5|(1 + 2)|3) 2 ^ 2mh 

+ 5 5 — hs_i (S34,ti23 512,556) 

l .[23](56)(4|(l + 2)|3) 3 (4|(2 + 3)|l) [23] (56) (4|(1 + 2)|3) 3 J V ; 



(24W4|f2 + 3)|l) -H-*^' -t^) 



[56] t 23i (24) (4|(2 + 3)|1) 

_ 1 [14] ({2| (1 + 4)(2 + 3)|5) 2 - (25) 2 g 23 ) 3 

2 (14)[23](56)(2|(1 + 4)|3)(2|(3 + 4)|1) 3 ^ 14 ' S23 ' S56 ^ 

1 [24] (5| (2 + 3)|4) 2 (2 S34 g 5 6 + ^12^234) (23)[16](3|(2 + 4)|1)(5|(2 + 3)|4) 



[23] [34] (56) ^34(11(3 + 4)12) i| 34 (4|(2 + 3)|1) 

+ 1 (23) [16] (3|(1 + 4)|6)(3|(2 + 4)|1) + 1 [14] (25) [16] (3|2|4) + ^ (23) [24] 2 (5|1 



2 <34> [56] *234 <4|(2 + 3)|1> 2 [34] i 234 <4|(2 + 3)|1) [23] t 234 <1|(3 + 4)|2) 

(3|2|1)(5|(2+3)|1) [46] 
[23] t 123 (4|(2 + 3)|1) 

I 1 r [14](3|l|6)(5|2|4) [14](5|4|2)(5|(2 + 3)|4)g 12 

2 (1|(3 + 4)|2)(4|(2 + 3)|1)L [34] + [23] [34] (56) 

^ ^ ((3|2|6)(3|(1 + 4)|2) + (3|1|2)(3|5|6)) + 2 ^ + 3 ^ 



(34)[56]^ 1 1 /x lv [23] 

: 56 + t 2 3< 

(4|123|4) 



+ ||||||||«5|2|4>( S56 -t 123 ) + (5|3|4)( S56 +t 234 )) + [14] (35) (2(3|4|6) + (3|2|6)) 



1 1 

+ 



2 (1|(3 + 4)|2)<4|(1 + 2)|3><4|(2 + 3)|1) L [23] ( < 5|(1 " ^ [W] ~ 2 [U] ^ 
+ [1|(3 + 4)2|4]«35) (4|1|6) + (45) <3|2|6>) + (5|3|6) (<3|2|1> s 13 + (3|4|1) s 12 ) 

+ [1|3(1 - 2)|6](3|(1 + 2)4|5) + (5|(3 + 4)|l)(3|4|6)s 2 3 + <5|4|1><3|(1 + 2)|6) s 13 
+ (3(3|2|l) + (3|4|l))(5|2|6) Sl4 -(5|3|l)(l|(3 + 4)|6)(3|(2-4)|l) 
+ 2 (3|4|1) [16] ((51) (a 12 + t 123 ) - (5|23|1)) - (3|2|1)(5|4|6) s 34 + 4 (5|123|1)(3|4|6) 
- (23) (4|(2 + 3)|1) (2 (5|(2 - 3)|4) [26] - (5|1|6) [24]) | I 3 3m (s 12 , s 34 , s 56 ) 

+ 2 (23) [24] (5|(2 + 3)|4)(5|(3 + 4)|2) L o (^f) _ ^ [14] [24] (15) (5|(2 + 3)|4) L 



[23] (56) *234 <1|(3 + 4)|2) s 34 [23] [34] (56) (1|(3 + 4)|2) s 56 

+ Aip 2 - 

The contributions where a scalar replaces the gluon in the loop are 



(1 



2e \-s 56 J 2' v 
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M 2 (l, 2,3,4) +M 3 (1, 2,3,4) + - 



1 (45)(3|(1 + 2)|4)(5|(2 + 3)|4) 



2 s 34 (56) (1|(3 + 4)|2)(4|(1 + 2)|3) 



In 



-S56 



1 (35) ((23) (45) - (25) (34)) 1 [24] (35) /(35) (5|(2 + 3)|4) \ 

2 (14)(34)(56)(4|(1 + 2)|3) + 2 [23](56)(1|(3 + 4)|2) V(34) + t 234 7 ^ ' 



(10.11) 

where M 2 and M 3 are given in eqs. (7.3) and (7.5). 

10.3 The Helicity Configuration q + q~ g + g~ 

Here we present the primitive amplitude A 6 (l+ , 2^~, 3 + , 4~). This amplitude has the same flip 
symmetry as the last one, as given in eq. (6.8). The tree amplitude is 



a tree 
^6 



— I 



[13] 2 (52)(4|(1 + 3)|6) (24) 2 [61](5|(2 + 4)|3) , [13] (24) [61] (52) 



[41] S 3i t 13i S 56 (23) S 34 t 234 S 56 

An alternate form with manifest k 5 \\ k e behavior is 



+ 



(23) [41] s 34 s 56 



A tree 
^6 



[13] 2 (52) (5|(2 + 4)|3) (24) 2 [61] (4|(1 + 3)|6) + [13] (24) [61] (52) 



L [34] [41] (56) i 134 i 234 (23) (34) [56] t 13A t 23A 
The results for the cut-constructible pieces are 



[41] (23) i 134 i 234 J 



(10.12) 



(10.13) 



V c 



1 



I 1 



(12)^[46] 2 tf 23 



(1|(2 + 3)|6) 2 (2|(1 + 3)|4) 2 



(23) [56] (1|(2 + 3)|4) 3 (3|(1 + 2)|4) (23) [56] (1|(2 + 3)|4) 3 (3|(1 + 2)|4) 



(10.14) 



Ls_™ (si4, t 423 ; S 23 , S56) 



+ 



+ 



(2|(1 + 3)|4) 2 (3|(1 + 2)|6)^ 



(23) [46] 2 ^ 23 

3 rrr^l /l l/o , o\M\ /0I/1 , oM A\ 3 



(23) [56] (1|(2 + 3)|4)(3|(1 + 2)|4) J [56] (1|(2 + 3)|4)(3|(1 + 2)|4) 



^ Ls 2 _ I 1 1 ' l (s 34 ,ti 23 ;si 2 ,S5 6 ) 



[23] 2 (15) 2 i 234 (1|(2 + 4)|3) 



(1|(2 + 4)|3) 3 (5|(3 + 4)|2)^ 



[34] (56) (1|(2 + 3)|4)(1|(3 + 4)|2) J [34] (56) t 234 (1|(2 + 3)|4)(1|(3 + 4)|2)' 



(24) 3 [16]' 



(23) (34) [56] t 234 (2|(3 + 4)|1) 



\ Ls 2 _™ h (si 2 ,t 234 ;s 34 ,s 56 ) 



+ 



+ 



(12) 2 (3|(1 + 2)|6)^ 



(23) (1|(2 + 3)|6)' 



(23) [56] (13) 2 (1|(2 + 3)|4)(3|(1 + 2)|4) [56] (13) 2 (1|(2 + 3)|4)(3|(1 + 2)|4) 

(5|(2 + 3)|4) 2 [23] 2 (5|(3 + 4)|2) 2 [34] 2 

(56) t 234 [24] 3 (l|(2 + 3)|4) (56) t 234 [24] 3 (1| (2 + 3)|4) 



Ls_x f^,^2i) 

\ — *123 — *123 / 



-S 2 3 



£34. 



■t234 ' — *234 



) 



+ TI$ m (s 12 , s 34 , s 56 ) + T/ 3 dm (s 14 , s 23 , s 56 ) + 2 



(5|1|3)(1|(2+4)|3)(5|(2 + 4)|3) 
[34](56)(1|(2 + 3)|4)(1|(3 + 4)|2) s 56 



+ 2 



(4|2|3)(5|(2 + 3)|4)(5|(2+4)|3) L o (^f 



(56) i 234 [24] (1|(2 + 3)|4) 
+ Mi(l,4,3,2) + flip 2 . 



S23 



+ 2 



(4|2|3)(5|(2 + 4)|3)(5|(3 + 4)|2) L o ' 



S34 



(56) t 234 [24] (1|(3 + 4)|2) 



•S34 



(10.15) 
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The coefficients of the triangle integrals are 



<24> 2 [16] (5|(2 + 4)|3) 1 [13] (24) (25) (5|(2 + 4)|3) 1 (24) 2 [16] (4|(1 + 3)|6) 

2 (23) [34] (56) t 234 2 (23) (34) [56] i 234 



(23) t 2 



234 



(45)(2|1|3)(2|(1 + 3)|6) 1 1 

(23) t 123 (1|(2 + 3)|4) + 2 (1|(2 + 3)|4)(1|(3 + 4)|2) 



(5|1|6)(4|2|3)(1|(2+4)|3) 

*234 



, (1|(2 + 4)|3)(5|(2 + 4) 1 3) 2 (2s 34 s 56 + S 12 t 234 ) , <5 1 1 1 3> 2 ( < 1 2> S 12 + (1|(2 + 3)4|2» 



[34] (56) il 34 



(23) [34] (56) 



+ (24) j [26] 2 (l|(2 + 3)|4) _ (24) (45) (1|(2 + 3)4(1 - 2)|6) + (12) [36 ] (5|(1 - 2)|3>(a 23 + *134) 



(23) (34) [56] 
(12) (45) 
(23) 



(23) (34) 
([36] (2s 24 + <5 34 ) - 2 [3|15|6]) 
1 



(23) [34] 



(1|4|6) «5|1|3> aw - (5|2|3) s 23 ) + 5 (1|241|3)(5|4|6) 

(2|14|3)(5|(1 + 3)|6) 



2(1|(2 + 3)|4)(1|(3 + 4)|2)(3|(1 + 2)|4)L 

+ (5|4|3)(1|(2 + 3)|6)(si3 - s 24 ) + (1|2|3) ((5|2|6) s 34 + (5|4|6) s 13 + 2 
+ 3 ((5|431|6) - (5|3|6) s 24 ) - 2 ((5|l|6)(s 14 + s 24 ) + (5|143|6) 



(23) 



(10.16) 



(1|(2 + 3)|4)(3|(1+4)|2) 
(4|(1 + 2)|3) 



+ 



2)|3)(5|(1+4)|6) + (5|1|3)(4|2|6) 

[13] (45) (2|(1 + 3)|6) 



(3|(1 + 2)|4) 



((5|2|4)(3|1|6)-(3|1|4)(5|(2-3)|6)) 



(10.17) 



<3|(l + 2)|4)t 123 



The contributions where a scalar replaces the gluon in the loop are 



1 



V sc = — 



2e \-s 56 
56 



1 

+ 2' 
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[24] (3|(2 + 4)|1) _ [23] 2 (24) 2 (5|(2 + 3)|4) 2 ^(f^f^) 

[14](3|(1 + 2)|4) 2 ^' ' ' ^ (56)[24]i 234 (l|(2 + 3)|4) t 2 34 

i (23) 2 [46] 2 tf 23 2m/t 

+ (23)[56] (1|(2 + 3)|4)(3|(1 + 2)|4)3 Ls -i ( Wm, s 12 , s 56 ) 

[23] 2 (15) 2 t 23 4 (l|(2 + 4)|3) 2mh , , 



[34] (56) (1|(2 + 3)|4)(1|(3 + 4)|2)3 



(23) (1|(2 + 3)|6)" / rr-" 1 , > , r p'«" , 

+ (23 ) [56](13 ) 2 (1|(2 + 3 ) |4 ) (3|(1 + 2 ) |4 ) LS - 1 ^'^ +T/3 ^'^4, , 56 ) + T/ 3 (-u,-*,,-*) 

+ [13] f(12)(3|(l + 2)|6) 2 L o(^f) (23)(1|(2 + 3)|6) 2 L o(^f) 



[56] <13> <3|(1 + 2)|4> V <3|(1 + 2)|4) s 12 <1|(2 + 3)|4> s 23 



1 



[23](24) 2 /[23](5|(2 + 3)|4) 2 L i(^) [34] (5| (3 + 4)|2) 2 L i (^f ) 



+ 



2 (56) [24] t 234 V <1|(2 + 3)|4) s 2 23 <1|(3 + 4)|2> sj 4 

(24) [46] 2 tj 23 L o (^f ) [23] 2 (24) (15) 2 t 234 L ° (^) 



[56] (1|(2 + 3)|4)(3|(1 + 2)|4) 2 s 56 (56) (1|(2 + 3)|4)(1|(3 + 4)|2) 2 s 3i 

+ 



L^ W J \ IV 1 /I / \ I V 1 J I / ^00 \ IV 1 /I / \ I V 1 J\ I ^oi 

1 [13] [23] (15) 2 t 23 4 <l|(2 + 4)|3) Z[13] L i(^) 2 L o(^) \ 

2 [34](56)(1|(2 + 3)|4)(1|(3 + 4)|2) V[23] i 2 34 (1|(3 + 4)|2) t 234 J 



+ 



I [13] [46] / ( 3 |(l + 4)|6)t 123 -(3|(l + 2)|6)t 134 X 

2 a|(2 + 3)|4)(3|(l +4)l2K3lfl + 2)14) V M ; ^ [141 [561 / V" S W 

23] (2|(3 + 4)|1) 
(3|(1 + 2)|4)A§ 



3 (( 45 ) ((2|3|6)5 23 - (2|5|6)5 56 ) - (14) [16] ((25) ~5 U - 2<2|36|5») 



1 1 



2 (3|(1 + 2)|4)A 3 
[36] (2|(3 + 4)|1) 



(25) (45) [16] (2|(l + 4)|6) x , ~ , 
-756^ + [14lT56l ) ([13] 523 " [1|(5 + 6)2|3] ) 



[14] [56] 

-M 2o (l,2,4,3)-M 3o (l,2,4,3) 



(56) [14] [56] 

([16]5 23 -2[1|45|6]) - 4<4|1|3) [16] (25) 



In 



1 (35)(4|(1 + 2)|3)(5|(2 + 4)|3) + 1 [46] (4|(1 + 2)|3) 



2 S34 (56)(l|(3 + 4)|2)(3|(1 + 2)|4) V"<W 2 [56] (3|(1 + 2)|4)(1|(2 + 3)|4) A 3 
/ (3|(1 + 2)|6)((24) <5 56 - 2 (2|13|4)) (12) [46] ([13] £ 56 - 2 [1|24|3]K 
V (34) [34] / 

+ \ mwrmmhw^ ( ((5|2|1) * 12 - (5|6|1) ^ w - [23] « 25) - 2 (2|1(3 + 4)|5 ») 

1 (2|(3 + 4)|1) / 1 4e (25) [16] £ 14 (45) (25) 
+ 2 (3|(1 + 2)|4)A 3 V2 [14] (23) + (23) (56) ^ M 1 

+ I 1 f (12) [13] 2 [46] 2 _ (25) (45) (2|(1 + 3)|4) + (24) (45) [46] ^ 



2 (3|(1 + 2)|4)(1|(2 + 3)|4) V [14] [34] [56] (23) (56) (34) J 

(24) 2 (5|(3 + 4)|2)(5|(2 + 3)|4) 1 (15) (45) (4|(1 + 2)|3) 

+ n ,ia\ lrn\ HUn , .Mol^l/o I oM 7\ + fll P2 • 



2 (23) (34) (56) i 234 [24] (1| (2 + 3)|4) 2 (34) (56) (1|(3 + 4)|2)(1|(2 + 3)|4) 

(10.19) 
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T = 3 
+ 



The three-mass triangle coefficients are 

gig [46] (56)<4|(1 + 2)|3) 
(3|(1 + 2)|4)(1|(2 + 3)|4)A| 
(12) [46] 



(12) [56] ((5|6|1) <5 56 - (5|2|1) 5 12 ) + t 124 ((2|5|6) 5 56 - (2|1|6) S 12 ) 

(3|(1 + 2 )|4><1|( 2 + 3)|4) A 3 [" [U] (56) ( 3 (4|(1 + 2)|3)(2|(1 + 3)|6) " (24) [36] { ^ ~ ^ 
+ (4|(1 + 2)|3)«5|2|1) 6 12 - (5|6|1) 5 5e + t 123 (5|(2 + 4)|1) - t 124 (5|(2 + 3)|1>) 

(4|(1 + 2)3|5) ([1|(3 + 4)(1 + 2)|4] - [1|23|4]) " 



+ (45)t 123 ([l|(3 + 4)(l + 2)|3]-[l|24|3]) - 

(2|1|3)(5|4|6) 
(3|(1 + 2)|4)(1|(2 + 3)|4) 



(3|(1 + 2)|4) 



+ 3 



+ 



(1|5|6)(4|(1+2)|3) 
(1|(2 + 3)|4)A| 
(15) 



(2|1|3)((5|2|1) 6 12 - <5|6|1><J 56 ) + <1|2|3>«5|6|2) <5 56 - (5|1|2) S 12 ) 



(2|(3 + 4)|1) - 
(1|(3 + 4)|2)_ 



(1|(2 + 3)|4)A 3 
<2|(3 + 4)|1> 



+ [23] 



[13] ((2|1|3)(4|5|6) + (4|2|3)(2|5|6) - (4|(1 + 2)|3)(2|(1 - 5)|6) 

[23] (14) <l|5|6)(* 23 4 -tl34) 



T = 3 



(1|(3 + 4)|2) 

(14) [23] (2|(3 + 4)|1) 
(3|(1 + 2)|4)A| 
(2|(3 + 4)|1) 
(3|(1 + 2)|4)A 3 



3(1|5|6)(4|2|3)-(1|2|6)(4|(1 + 2)|3) + 



(1|(3 + 4)|2) 



(S 14 (2|5|6) - s 56 (2|3|6)) (5|4|1) - - $ 56 <2|5|6><5|6|1> 



-(4|5|6)((5|6|3)+2(5|2|3))-(5|2|3)(4|1|6) 



(10.20) 
(10.21) 



11. Contributions with the Vector Boson coupled to a Fermion Loop 

We now consider the remaining contributions where the vector boson (7* , Z) couples directly 
to a quark loop, as depicted in figs, le and f. The contributions proportional to the vector and axial- 
vector couplings of the quark to the vector boson are separately gauge invariant and have different 
symmetry properties, so we separate the two contributions. Both contributions are infrared and 
ultraviolet finite, because there is no tree-level coupling between the vector boson and any number 
of gluons. Therefore for each amplitude in this section V = 0, and we just give the finite (F) terms 
in eq. (6.1). 

For the vector coupling case, Furry's theorem (charge conjugation) implies that only box 
diagrams contribute. The three box diagrams are shown in fig. le. A simple way to understand 
the cancellation of triangle diagrams, depicted in fig. 9, is that under reversal of the fermion loop 
arrow the sign of the diagrams flip so that there is a pairwise cancellation. (This argument also 
relies on the existence of only one £[/(iV c )-invariant combination of two gluons, that proportional 
to 5 ab .) For the axial- vector case, Furry's theorem does not apply, and triangle diagrams such as 
those in fig. 9 do contribute. 
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Figure 9. Two examples of triangle diagrams that contribute to the axial-vector coupling 
case, but cancel in the vector coupling case. 



For a massless isodoublet the axial-vector contribution exactly cancels. Since we take the 
u,d,s,c,b quarks to be massless, but not the t quark, only the t,b pair survives the isodoublet 
cancellation in the loop. As a uniform approximation we keep all terms of order 1/to 2 but drop 
terms of order \jm\ and beyond. For the vector coupling, the contribution of the top quark loop 
decouples rapidly, leaving only terms of order \jm\ which we drop. For the axial- vector part, 
terms of order 1/m 2 appear which we keep. The 1/m 2 terms are rather easily obtained (especially 
when compared to a diagrammatic calculation) from an effective Lagrangian (or operator product) 
analysis, where the coefficients are fixed via the collinear limits. 



11.1 The Helicity Configuration q + q g + g 



We now present the results for the primitive amplitudes Aq s (1+ , 2^ , 3 + , 4 + ), A§ (1+, 2^ , 3 + , 4 + ), 
^r(l Q + ,2r,3+,4+),and ^ Bl (l+ 2r, 3+, 4+). 

The various contributions are 
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(11.1) 



«12 
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(11.2) 
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™*_ 1 ((23) (45) + (24) (35)) (25) 2me , 
2 (12) (56) (34) 3 V 
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,S 5 <i) + 
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S34 
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(25) 
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134 
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(24) 



(11.3) 



^7>ax,sl 
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(13) (23) s 56 
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tl23 \ 
S56 / 



S56 



+ 



24 ml j 



+ exch 34 . 



(11.4) 



11.2 The Helicity Configurations q + q g + g and g + g g g^ 



It is convenient to quote the helicity amplitudes yl6 ,ax (l+, 2^, 3+, 4") and Al ,ax (l+ , 2^, 3", 4+) 
together because they are rather similar. In fact for the vector case Furry's theorem (charge 
conjugation) implies that the two A$ amplitudes are identical up to the relabelling 3^4, 



^(l+,2r,3-,4+) = ^(l+,2r,3+,4-)| 



3^4 



(11.5) 
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For A%(1+, 2^ , 3+ , 4") we have 
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(11.6) 
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(11.7) 



The full axial-vector amplitude does not obey a relation as simple as eq. (11.5) under 
exchange of the two gluons. However, the more complicated parts of 4| x (the box diagram contri- 
butions) do obey such a relation, with an additional minus sign due to the 75 insertion. Thus it 
is convenient to separate the vlg x contributions into a common part (7 ax obeying the relation, plus 
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additional terms with no special symmetry. The common part is 
1 (2|(1 + 3)|4) 2 (3|(1 + 2)|6) 2 - <23> 2 [46] 2 tj 23 2m , 



2 (12) [56] (3|(1 + 2)|4)4 
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-6 



(3|(1 + 2)|4) 2 A 3 2 i 123 (3|(l + 2)|4) 

12] (2|(1 + 3)|6)«25) 5 M - 2<2|16|5» (4|(1 + 2)|3) [13] [46] (2|(1 + 3)|6) 
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(2|(1 + 3)|6) 2 ,x 
(12)[56](3|(1 + 2)|4) 2 V-W 
(24) [36] / (2|4|6)^ 34 _ (24) (35) <5 56 _ [13] [46] 5 12 _ <5|3|1) (24) (35) \ 
+ (3|(1 + 2)|4) V(12)[56] A 3 (12) (34) A 3 [34] [56] A 3 A 3 ^ (12) (34) s 56 J ' 

(11.8) 

where the operation exchi 6j25 is to be applied to all preceding terms within the brackets. 
For ^ x (l+,2^,3+,4-) we have 

(24)(l|4|6)(2|(l+3)|6)t 123 L i(^S) (2|(1 + 3)|6)(3|(1 + 2)|6) [13] L o(^f 
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+ (13)(34) S56 (3|(1 + 2)|4) 12m 2 (13) (34) s 56 + P2 ' 



(11.9) 



where the flip 2 operation is to be applied to all terms including the C ax terms. 
Similarly, for ,4g x (l+, 2^, 3~, 4+) we have 

[14] 2 (45)(5|(2 + 3)|l)t 123 L i(^S) (5|(2 + 3)|1)(5|(1 + 2)|3)(23) L o(^) 



-C^L .„ + 



+ 
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[14](4|(2 + 3)|1)(5|(2 + 3)|1)(5|(1 + 2)|3) 
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(11.10) 

The subleading-color axial-vector contributions do obey a 3 — 4 symmetry relation, 

^r Sl (1^2r,3-,4 + ) = Ar sI (l Q + ,2r, 3 + ,4-)| 3 ^ 4 . (11.11) 
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We give here ^ x ' sl (l+, 2r, 3+, 4"): 



, sl _ (24) (45) <2|(1 + 3)|6) 
(13) (23) S56 



S56 



+ 



24m t 2 j 



+ flip 



2 ' 



(11.12) 



12. Simplified Versions of Four Quark Amplitudes 

In a previous paper we have presented the one-loop helicity amplitudes for e + e~ — > qqQQ 
[10]. Here we present simplified versions of these amplitudes. These versions have fewer spurious 
singularities and are therefore (slightly) better for implementing in a jet program. We have verified 
that the two forms are numerically identical. (In contrast to the e + e~ — > qqgg case, here there is 
not much to gain from decomposing the amplitude into a scalar contribution and a cut-constructible 
part.) 

The primitive amplitudes and A s & are proportional to tree amplitudes and are given by 



^■ +± (1, 2, 3, 4) = c r 4 ree (l+ 2±, 3J, 4,; 5 e " , 6+) 



3e V-S23 
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(12.1) 



where ^4g ree is given below. 

The top quark vacuum polarization contribution is 



4< +± (l, 2, 3, 4) = -1 c r § j4|-(1+ 2± 3* 4,; 5r , 6+) , 



15 ml 



(12.2) 



neglecting corrections of order l/m| and higher. The remaining contributions are decomposed 
further into divergent (V) and finite (F) pieces according to eq. (6.1). 

12.1 The Helicity Configuration q + Q + Q~q~ 

We first give the primitive amplitude A§ + (1, 2, 3, 4). This amplitude is odd under the operation 
flipj^ defined in eq. (6.7). The tree amplitude for this helicity configuration is 



< ee < ++ (l,2,3,4)=* 
The loop amplitude is 



[12] (54) (3|(1 + 2)|6) + (34) [61] (5|(3 + 4)|2) 



S23S56*123 



523^56^234 



(12.3) 



63 
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(12.5) 

where is to be applied to all the preceding terms in F ++ . The structure of this amplitude is 
already rather simple, and it is unchanged from ref. [10]. 



12.2 The Helicity Configuration q + Q~Q + q 



We now give the result for A§ (1, 2, 3, 4). This amplitude is odd under the same flip symme- 
try (6.7) as Aq + . The tree amplitude is 
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Note that A^ ee ' + (1, 2, 3, 4) = -A" ee ' 3, 2, 4). The simpler form for the loop amplitude is 



i tree, H — h 



^+-(1,2,3,4) = y++(l,2,3,4), 

64 



(12.7) 
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2 S i2(-(25)[16] + (15)[26]|||^)]ln(5 S ) 

- fliPi , 

where flip 1 is to be applied to all the preceding terms in F + ~ . 

12.3 Subleading-Color Primitive Amplitude 

The primitive amplitude Af(l, 2, 3, 4), contributes only at subleading order in N c . The "tree 
amplitude" appearing in eq. (6.1) is 
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and satisfies A^ ree ' sl (2,3, 1,4) = -A 



tree, ++ / 



6 y^^^^j — (1)2,3,4). 

To describe the loop amplitude, we use the operations flip 3 and exchi 6j25 defined in eqs. (6.9) 
and (6.11). The loop amplitude is given by 
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where exchie j2 5 is to be applied to all the preceding terms in F sl , but flip 3 is to be applied only 
to the terms inside the brackets ([ ]) in which it appears. The three-mass triangle coefficient T is 
given by 

l (3s 3 4<534-A 3 )(t 123 -t 124 )(2|(3 + 4)|l)(5|(3 + 4)|6) 1 s 34 (t 123 - t 124 ) (25) [16] 
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12.4 Axial- Vector Quark Loop 

The axial- vector quark-triangle contribution ^4| x (l, 2, 3, 4) is easily obtained from the fully 
off-shell Zgg vertex presented in ref. [49]. The infrared- and ultraviolet-finite result is 
, axn 9 „ _ 2i f(m t ; s 12 ,s 34 , s 56 ) - f(m b ; s 12 ,s M , s 56 ) ( [63] (42) (25) [61] [13] (45) 
6( '' 3 ' } " Wf ^ V (12) [12] 

+ (1 ~ 3, 2^4), 

(12.13) 

where the integral f(m) is defined in eq. (11.14). We need only the large mass expansion (for 
m = m t ) and the m = limit (for m = m^) of this integral; these are presented in eq. (11.15). 

13. Summary and Conclusions 

In this paper we presented explicit formulas for all one-loop helicity amplitudes which enter into 



66 



numerical programs for the next-to-leading order QCD corrections to e + e~ — > (j*,Z) — > 4 jets. 
These amplitudes have already been incorporated into the first such program to be constructed [12]; 
the result has been a significant reduction in theoretical uncertainties in the four-jet cross-section 
and associated quantities. With appropriate modifications to the coupling constants, these same 
amplitudes enter into the computation of next-to-leading order contributions to the production of 
a vector boson (W, Z, or Drell-Yan pair) in association with two jets at hadron colliders, and 
three-jet production at ep colliders. 

Following the methods reviewed in ref. [25], we have obtained the loop amplitudes by demand- 
ing that their functional forms satisfy unitarity and factorization. This approach makes use of a 
color decomposition [32] into gauge invariant primitive amplitudes [19,10] which are expressed in 
terms of a helicity basis [31]. The color decomposition limits the analytic functions that may ap- 
pear, greatly simplifying the reconstruction of the amplitudes from their analytic properties. The 
helicity basis results in relatively compact expressions, from which spurious poles can be systemat- 
ically removed to further simplify the results. A further advantage of the helicity basis is that one 
retains all spin information. 

As a check on the methods we have verified that the amplitudes are numerically identical 
to ones obtained by a (numerical) Feynman diagram calculation that we have performed. This 
diagrammatic calculation made use of a number of string-motivated ideas reviewed in ref. [25]. We 
have also numerically verified that the typed form of the amplitudes appearing in this paper agree 
with our Maple and Mathematica expressions for the same quantities. 

We expect that the amplitudes presented in this paper should lead to an improved knowledge 
of QCD predictions for a wider class of observables, and thus of the QCD background to searches 
for new physics in various processes. 
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Appendix I. One-Loop Integrals 



LI General Properties and Reduction Formulae 

The loop momentum integrals that appear in either a Feynman diagram or cut-based analysis 
are of the form 

/AD p pan . . . pa m 
7^77 i , (I-l) 

where the Ki are external momenta, or sums of external momenta. The momentum £ flows through 
the propagator between external legs n and 1. For convenience we define 

I°[£ ai ---£ a ™\ =i(-l) n+1 (4vr) D / 2 jf [£ ai ■ ■ ■ £ a ™} . (1.2) 

Integrals with powers of loop momenta £ a in the numerator are known as tensor integrals; integrals 
with no powers of loop momenta in the numerators are known as scalar integrals and are denoted 
by I® = When the superscript D is omitted below, D is implicitly equal to 4 — 2e. 

In general, any one-loop amplitude may be expressed as a linear combination of bubble, triangle 
and box (i.e. two-, three- and four-point) scalar integrals. This follows from the Passarino-Veltman 
reduction [39] of any tensor n-point amplitude with m powers of loop momenta to a linear com- 
bination of n — 1 and n-point integrals with no more than m — 1 powers of loop momenta. The 
resulting scalar integrals with n > 4 legs can be further reduced to scalar box integrals using an 
additional set of reduction formula? [40,41,42]. 

It is useful to review briefly the conventional reduction of integrals to explain the origin and 
types of determinantal poles that can appear. Consider the five-point tensor integral with the 
kinematic configuration depicted in fig. 6. Following the Passarino-Veltman reduction technique, 
we expand the loop momentum in terms of four independent external momenta, 

4 

h[£ a H a * ...£ a >\ =J>r4fc (1-3) 
k=i 

where 

Pi=Ki, p 2 = K 1 + K 2 , p 3 = K 1 + K 2 + K 3 , p 4 = -K 5 , (1.4) 

and we have suppressed the Lorentz indices a 2 ■ ■ ■ ctj on the right-hand side of eq. (1.3) in A^. The 
functions Ai are found by first contracting eq. (1.3) with the momentum sums pi, generating the 
four linearly independent equations, 

4 

2I 5 [£ ■ Pl £<**... r*\ =J2*ikA k , (i = l,2,3,4), (1.5) 

fe=i 
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where tik = 2pi ■ pj~. Using 2£ ■ pi = — (£ — pi) 2 + £ 2 + p 2 , and recognizing that (£ — Pi) 2 and £ 2 
are inverse propagators of the pentagon integral, one can reduce the left-hand-side to a sum of two 
four-point integrals plus a five-point integral with one less power of loop momentum. Inverting Uk 
we have 

4 

Ai = Y J t 7k ( J i fe +1) [i a2 i as ■■■£ aj ] -l1\l a n a3 ■■■£ a i] +p\h[£ a2 £ a2 ■■■^]) , (1.6) 
fe=i 

where /I*' 1 denotes the box integral that is obtained from the pentagon integral by canceling the 
propagator factor l/£ 2 . The coefficients contain in their denominators the pentagon Gram 
determinant A 5 given in eq. (1. 10). 

The reduction of tensor boxes, triangles, and bubbles is similar except that one must also 
include the Kronecker-delta in the expansion of the integral function since there are less than four 
independent momenta. More powerful reduction techniques which lead to square-roots of Gram 
determinants in denominators, instead of single powers, have also been developed [40]. 

The scalar pentagon integrals may be reduced to scalar box integrals using the scalar integral 
recursion formula [40,41,42] (valid for n < 6) 



i = l 



(1.7) 



where 

j = l i=l i,j = l 

The matrix Sij is defined in eq. (4.4) and the integral I® =6 ~ 2e is the scalar n-point integral evaluated 
in (6 — 2e) dimensions. Observe that for n = 5 the prefactor of j^= 6 - 2e [ s Q f 0(e); since i® =Q ~ 2e 
is finite as e — ► we may drop this term. This equation is also useful for rewriting D = 4 — 2e 
box integrals as D = 6 — 2e box integrals, which turns out to be a convenient way to represent 
the amplitudes. Indeed, the Ls_i functions defined in appendix II are D = 6 — 2e box integrals 
from which a simple overall kinematic factor has been removed. The explicit forms of the higher- 
dimension box integrals may be obtained by solving eq. (1.7) for if =6_2e 

For the case of a single power of loop momentum we may avoid a Gram determinant in the 
denominator by making use of the reduction formula, 



1 n 



2 

i=2 



.3 = 1 



(1.9) 



where pi = K\ + K 2 H h -fQ and the Kj are the external momenta of the integrals. For pentagon 

integrals (n = 5) the i^= 6 ~ 2e term may be dropped since it is of 0(e). 
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Since the integral recursion formula (1.7) contains inverses of the matrix S defined in eq. (4.4), 
the amplitudes will contain poles in det S. As discussed in the text, these poles do not correspond 
to the propagation of physical states, and so their residues in a full amplitude must vanish. 

In summary, when reducing tensor and scalar integrals to linear combinations of scalar box, 
triangle and bubble integrals, one encounters a set of spurious determinantal poles. Some of these 
spurious poles, such as the pentagon Gram determinant (see section 4.2), are artifacts of the 
reduction procedure, but others are inherently part of the amplitude when it is expressed in terms 
of logarithms and dilogarithms. 

1.2 Determinants Appearing in Amplitudes 

Here we explicitly list the determinants that appear in the denominators of coefficients in the 
general integral reduction formulae of the preceding subsection. As discussed in section 4, knowledge 
of these determinants — and how they can be factored — is useful in deciding which spinor factors 
to multiply and divide by in order to simplify the integral reduction of specific cuts. 

The Gram determinants are defined in eq. (4.2). The explicit forms of the Gram determinants 
associated with pentagon and box integrals whose external legs follow the ordering 123456 is 

A, — _ s 2 „2 _ 2 f 2 _ f 2 f 2 _ 2 f 2 _ 2 2 _ 2 2 
^5— *23 6 34 6 34 L 123 ''123 ''234 6 12 ''234 6 23 6 56 6 12 6 23 

2 2 2 2 2 o 2 

+ 2 S 34 £i23 S23 + 2 t 12 3 *234 S34 + 2 S 12 t 2 34 ^23 + 2 S56 S 23 S34 + 2 S56 S 23 S 12 + 2 S 23 «12 «34 

— 2 ti23 *234 S23 S34 ~~ 2 t 2 34 «12 323 S34 ~~ 2 t 234 S34 t i2 3 S 12 — 2 t i23 «12 «23 S34 + 2 t 234 *123 s 12 

— 2 ti23 ^234 Sl2 S23 — 2 S56 «23 S34 ^123 + 2 S56 «23 *123 *234 ~ 2 S56 «23 S 12 ^234 + 4 S56 S 2 3 Si 2 «34 , 

(1.10) 

-2S23S34S24 , 



A 



(i 



A 



A 
A 



(3 



(4 



A 



(5 



(1.11) 



2 S34(S12 S56 — ^123 il24) j 
2 «14(S23 £56 — *123 ^234) , 
2 S12(S34 S56 — tl34 ^234) j 
-2si 2 S 23 Si 3 , 

where the superscript on the box Gram determinants labels the propagator that has been canceled 
in the pentagon integral to obtain the box. Similarly, the triangle Gram determinants are labeled 
by a pair of indices, for the two canceled propagators, 



A (l,2) 
^3 


- s 2 

— *34 > 


A (1.3) 
^3 


— (^234 — 


« ^2 a (1,4) 
s 23j , ^3 


— (^234 — 


a ^2 A (1 ' 5 ) - c 2 
s 34j , A3 - S 23 , 


A (2,3) 
^3 


= (*123 


- S56) 2 , 


^3 


2,2,2 

— s 12 + s 34 -t- s 56 


2s i2 S34 - 


- 2s 34 s 56 - 2s 56 si 2 , 


a(2,5) 
^3 


= (*123 


- S 12) 2 , 


A (3,4) 
^3 


= (^234 — S 5& ) 2 , 


A (3,5) 
^3 


= (*123 — S23) 2 , 


a(4,5) 
^3 


- s 2 

— S 12 • 













(1.12) 
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The bubble Gram determinants are trivial. As discussed in section 4, the pentagon Gram determi- 
nant does not appear at all in the cut calculations (assuming the integral reductions are performed 
as discussed). However, the triangle and box Gram determinants do appear in the final results, 
when expressed in terms of logarithms and dilogarithms. 

Equations (1.7) and (1.8) show that another source of spurious singularities is the determinant 
of the matrix Sij. Again labeling the pentagon and box matrices by the canceled propagators we 
have 

det S = -^ s 12 S 2 3 S34(S23 «56 ~ *123 *234) (1-13) 

det = — s 2 23 s| 4 , det = ^ ^ ? det 5 (3) = _ tl23 t 234 ) 2 , 

io ID (1.14) 

det S^- ^ = — s 12 t 234 , det = — s 12 s 23 . 

Since the primitive amplitudes contain loop integrals where the external legs follow the order- 
ings 123564 and 125634 instead of 123456, an additional set of possible spurious poles are obtained 
from the above set via the relabelings of external legs: 

1234 ^4123, 1234 ^3412. (1.15) 

Observe the appearance of the combinations 

812 S56 " tl23 *124 = ~(4|(1 + 2)|3> <3 1 ( 1 + 2)|4) , 

S23 «56 " *123 *234 = -<1|(2 + 3)|4) <4|(2 + 3)|1) , (1.16) 
^34 S 56 " 1 134 *234 = - (1 1 (3 + 4) |2) <2| (3 + 4) 1 1) , 

which we have factored into products of spinor strings. The above factored form is quite useful in 
simplifying the cuts. As discussed in section 4.2, we insert such 'back-to-back' factors into tensor 
integrals by hand in order to help simplify the expression by forming appropriate spinor strings 
from which inverse propagators can be extracted. 

Although the three-mass triangle Gram determinant A 3 = Ag 2 ' 4 - 1 cannot be factored simply 
like eq. (1.16), there are many kinematic combinations that vanish whenever A 3 does. Using the 
fact that the three four- vectors k\ + k 2 , k 3 + /c 4 and k 5 + k 6 all become proportional in the limit 
A 3 — > 0, and relations like 5$$ = 2 (k\ + k 2 ) ■ (Aft + fc 4 ), it is easy to verify that the expressions 

(1|(3 + 4)|2), (14)«5 56 - 2 (1|23|4), (1|2|4) «5 12 - (1|3|4) S 3A , 

(1.17) 

(l|(l + 2)(3 + 4)|6)-(l|(3 + 4)(l+2)|6), 

plus their complex conjugates and a variety of permutations of them, all vanish in this limit. In 
the amplitudes we present, the appearance of such combinations in the numerator of coefficients, 
when A 3 appears in the denominator, alleviates the amplitudes' spurious singularities as A 3 — > 0. 
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Appendix II. Integral Functions Appearing in Amplitudes 

We collect here the integral functions appearing in the text, which contain all logarithms and 
dilogarithms present in the amplitudes. Most of the functions have already appeared in previous 
papers [17,10], but for completeness we list them all in this appendix. Except for the contribution 
of the top quark to vacuum polarization contributions and to the axial- vector contribution A^, all 
internal lines are taken to be massless. The following functions arise from box integrals with one 
external mass: 

L °W = T37' l() = l_ r ■ 

Ls_i(ri,r 2 ) = Li 2 (l - n) +Li 2 (l - r 2 ) + lnr 1 lnr 2 - — , 

1 6 (HI) 

Ls (ri,r 2 ) = Ls_i(ri,r 2 ) , 

(1 - n - r 2 ) 

Lsi(ri,r 2 ) = ^ _ — _ — ^ [Ls (ri,r 2 ) + L (r 1 ) +L (r 2 )] , 
where the dilogarithm is 

Li 2 (x) = - fdy^A. (II.2) 

Jo v 

The function Ls_i is simply related to the scalar box integral with one external mass, evaluated 
in six space-time dimensions where it is infrared- and ultraviolet-finite. The above functions have 
the property that they are finite as as their denominators vanish. Generalizations of the Lso and 
Lsi functions to the case of box integrals with two or more external masses have been presented in 
ref. [38]. 

The box function analogous to Ls_i, but for two adjacent external masses, is 
Ls 2 jf {a, t; ml m \) = -U 2 (l - ^] -U 2 ( I - ^ -\ In 2 In 



+ 



1 2 o m\m\ 
-{s-m{- m 2 ) H — 



^ — 1 1 2 \ —mi I V — mi 
/| m ( S ,m 2 ,m 2 ), 



(II.3) 

where l| m is the three-mass scalar triangle integral. This integral vanishes in the appropriate 
'back-to-back' kinematic limit. We also employ a version of this box function with /| m removed, 
' 2mh ' ■ 2 t: fi m l\ r; fi m 2 A 1 ^f~ s \ 1 



(II.4) 

The label l h' refers to the fact that this integral is relatively hard to obtain [42]. 
The box function with two non-adjacent external masses is 

2me/„ . 2 „9n _ T; f, ™<i\ j ■ {, mj\ ( m\\ T . (_ m\ 



Ls 2 _r(^t;m 2 ,m 2 ) = -Li 2 (l--ij - Li 2 \\ - -± j - Li 2 ^ 1 - ^ ) - Li 2 ( 1 - 

mimW 1 _ i ( ~s 
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(II.5) 



In this case the label 'e' refers to the fact that this integral is relatively easy to obtain. This integral 
vanishes as s + 1 — m\ — m| — > 0. 

The kinematic region in which an Ls_i function vanishes always turns out to be related to 
the spinor product (or string) required to reduce the corresponding box integral: (24) for 
(3|(1+2)|4) for lf \ (14) for lf \ (1|(3+4)|2) for and (13) for lf ] (or their complex conjugates). 
Curiously no single helicity amplitude contains both Ls 2 _ I 1 l/l and Ls 2 .™ 6 functions simultaneously. 

The analytic properties of these integrals are straightforward to obtain from the prescription 
of adding a small positive imaginary part to each invariant, — > s^- + is. One expands the 
logarithmic ratios, ln(r) = ln(£p-) = ln(— s) — ln(— s'), and then uses 

In(-s-ie) = ln|s| -otG(s). (II.6) 

where 6(s) is the step function: B(s > 0) = 1 and 6(s < 0) = 0. The imaginary part of the 
dilogarithm Li 2 (l — r) is given in terms of the logarithmic ratio, 

ImLi 2 (l — r) = — ln(l — r) Imln(r). (H-7) 

For r > the real part of Li 2 (l — r) is given directly by eq. (II. 2). For r < one may use [50] 

7T 2 

ReLi 2 (l-r) = — - In \r\ In |1 - r\ - ReLi 2 (r) , (II.8) 

with ReLi 2 (r) given by eq. (II. 2). 

The analytic structure of /| m is more complicated [51,52,42], and the numerical representation 
we use depends on the kinematics. The integral is defined by 

f 1 1 

^! m (si2,s 3 4,s 5 6) = / d 3 a i 5{l-a 1 -a 2 -a 3 ) . (II.9) 

Jo -si 2 aia 2 - s 34 a 2 a 3 - s 56 a 3 ai 

This integral is symmetric under any permutation of its three arguments, and acquires a minus 
sign when the signs of all three arguments are simultaneously reversed. Therefore we only have to 
consider two cases, 

1. The Euclidean region si 2 ,S34,S56 < 0, which is related by the sign flip to the pure Minkowski 
region (si 2 ,s 3 4,s 56 > 0) relevant for e + e~ annihilation. Here the imaginary part vanishes. This 
region has two sub-cases, depending on the sign of the Gram determinant A 3 (si 2 , s 34 , s 56 ) defined 
in eq. (6.6): 

la. A 3 < 0, 
lb. A 3 > 0. 

2. The mixed region s\ 2 , < 0, s 3 4 > 0, for which A 3 is always positive. 
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In region la one may use a symmetric representation found by Lu and Perez [51], which is 
closely related to that given in ref. [42]: 



r3m 



Cl 2 2tan" 



'12 



+ Cl 2 2 tan 



i 



J34 



+ Cl 2 2 tan 



J56 



where the 5ij are defined in eq. (6.6) and the Clausen function Cl 2 (x) is defined by 



C\ 2 (x) ee 



sin(nx) 



n=l 



dt ln(|2sin(t/2)|). 



(11.10) 



(11.11) 



In regions lb and 2 a convenient representation is given by Ussyukina and Davydychev [52], 



r3m 
J 3 



1 



Re 



/A3 
i7re(s 34 



2 (Li 2 (-px) + Li 2 (-py)) + ln(px) ln(py) + In (|) In ^ - 
(<5i2 + y/^)(5 56 +y^ 



1+py 



+ px / 3 



In 



(<5i2 - y/£d(he - V^s, 



where 



2s 



56 



Sl2 S34 

x = — , y = — , p = 

«56 s 56 (5 56 + vA 3 



(11.12) 



(11.13) 



Finally, in the top quark contribution to A§ x the combination f(m t ) — /(m-t,) appears, where 
f(m) is the integral 



/(m;si 2 ,s 3 4,s 56 ) 



f 

Jo 



d 3 ai 5(1 — a\ — a 2 — a 3 ) 



a 2 a 3 



(11.14) 



m' - s 12 aia 2 - s 3i a 2 a 3 - s 5 6a3«i 
This integral is complicated for arbitrary mass m; however, the large and small mass limits of it 
suffice for m t and mj respectively. For m = m t we simply Taylor expand the integrand in 1/m; for 
m = mi we set raj, to zero, and reduce /(0) to a linear combination of the massless scalar triangle 
integral I 3 m given above, logarithms and rational functions. We get 



/(w-t;si2,s 34 ,s 56 ) 
/(0;si 2 ,s 34 ,s 56 ) = 



i (2S 34 + 5i2 + S56) 

~r a t 



24m? 



360m- 1 



3s 3 4<^34 1 \ ,-;;„, , 

~T Si 2 S 56 / 3 (S12, S34, «56) + 



+ 



A 2 
^3 

3si 2 <5 



3S5fi(5; 



56"56 



A 2 
^3 



2A, 



S12 In 



12012 



A 2 
^3 



2A, 



s 56 In 



—■§56 
-S34 



"34 

2A 3 



Note that the limit where one of the invariants vanishes, 

f(m t ; 0, ti2 3 , s 56 ) - /(0; 0, t m , 550) = - 
appears in two amplitudes, eqs. (11.4) and (11.12). 
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(^) 



«56 



+ 



24 m 2 ' 



-■S12 
-S34 



(11.15) 



(11.16) 



Appendix III. Spinor Identities for Simplifying Spurious Poles 



In order to simplify or remove spurious poles a number of spinor identities are of great utility. 
In this appendix we collect these identities. As discussed in section 4 some of the spurious poles 
can be completely removed from the amplitudes, but others are an inherent part of the amplitude 
when it is expressed in terms of logarithms and dilogarithms. 

The most difficult spurious poles to simplify are those appearing as coefficients of three external 
mass triangle integrals and associated logarithms. In simplifying these spurious poles, it is useful to 
have identities for rewriting commonly occurring expressions to make one behavior or another more 
manifest. The following identities are used in different 'directions', depending on the situation. 

Particularly important spurious singularities are where A3, (3|(1 + 2)|4) and t\23 — S12 (and 
the label-permuted objects) vanish. In order to obtain relatively compact expressions it is essential 
to simplify these spurious poles. 

One important quantity is (ti23<^34 + 2si 2 s 56 ), which appears as a factor in the coefficient 
of /| m (si2, S34, s^q) within the function Ls^_™ h (s34, £123; S12, «56) defined in eq. (II. 3). If we shift 
from using the function Ls 2 _^ h in an amplitude, to Ls_ 1 instead (or vice- versa), then we also 
have to shift the coefficient of the three-mass triangle by a quantity which contains this factor, 
as well as the spurious singularity factor (3|(1 + 2)|4) typically associated with the coefficient of 
hs^ h (s 3i , i m ; S12, s 56 ). The identities 



are useful in manipulating leading 1/(3|(1 + 2)|4) singularities in order to remove, for example, 
'extra' l/ii.23 poles. Also note that 




(III.l) 



<$34 + *123 + tl24 — , 



$12 + ^234 + tl34 — 



$12 + #34 + 2s 56 = . 



(III.2) 



The identities 



Sl2 Sl2 + «34 £34 + s 56 8 56 = A 3 , 



(5 12 ) 2 = 4s 34 s 5 6 + A 3 , 

S56 #56 = \ $12 $31 + \ A; 



3 , 



(III.3) 
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are useful in simplifying the I/A3 poles. Also useful in this regard are the identities, 



S12 (35) [46] = (34) [46] 2 (65) + [43] (35) 2 [56] + (3|(1 + 2)|4)(5|(3 + 4)|6) , 

<5 34 (25) [16] = (21) [16] 2 (65) + [12] (25) 2 [56] - (2|(3 + 4)|1)(5|(3 + 4)|6) , (III.4) 

5 56 (32) [41] = (34) [41] 2 (12) + [43] (32) 2 [21] - (3|(1 + 2)|4)(2|(3 + 4)|1) . 

where we show a few permutations. A consequence of the last identity is 

(A 3 + 4s 12 s 56 ) (25) [16] = <5 34 ((21) [16] 2 (65) + [12] (25) 2 [56]) - <5 34 <2|(3 + 4)|1)(5|(3 + 4)|6) . 

(III.5) 

Another related identity is (see also eq. (5.21)) 

[64] (32) t 123 = (2|(1 + 3)|6)(3|(1 + 2)|4) + (3|(1 + 2)|6)(2|(5 + 6)|4) . (III.6) 



Appendix IV. e + e — > qqg Helicity Amplitudes 



In this appendix we present the e + e~ — ► qqg primitive amplitudes, which appear in the 
collinear limits discussed in section 3.2. These amplitudes were first calculated in the spinor helicity 
formalism by Giele and Glover [7]. We present these amplitudes in the same primitive amplitude 
format as for the e + e~ — > qqgg amplitudes discussed in section 2.3, including a separation into 
cut-constructible and scalar pieces. The parent diagrams associated with these amplitudes are 
depicted in fig. 10. Using parity and charge conjugation invariance, there is only one independent 
helicity configuration for each of two independent color configurations. 

These amplitudes appear in the singular collinear limits of the e + e~ — > qqgg and e + e~ — ► 
qqQQ amplitudes, except for the k$ \\ channel where the lepton pair becomes collinear. For 
the primitive amplitudes with closed fermion loops, the amplitudes e + e~ — > ggg also appear in 
the collinear limits. Although we do not present these amplitudes here, they may be obtained 
from ref. [53] after using the spinor helicity representation for the polarization vectors. One may 
also obtain the vector coupling results from the qqggg amplitudes of ref. [19] after summing over 
permutations of the quark lines, which effectively removes their color charge. 
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(a) A 5 (l q , 2, 3 fJ 



(6) A 5 (l q , 2 rJ , 3) 



(c)Al sJ (l q , 2 f7 , 3) 




(d) A?(l t , 2 q , 3 



Figure 10. Parent diagrams for the various two quark and one gluon primitive amplitudes. 
Straight lines represent fermions, curly lines gluons, and wavy lines a vector boson (7* or Z). 

First consider the primitive amplitude ^4s(l+, 2 + , 3^~, 4j , 5+), which contributes at leading 
order in 1/N C . The tree amplitude in this case is 



Af ee = -1 



(34)< 



(12) (23) (45) 

The results for the cut-constructible and scalar pieces are 



(IV.l) 



v c 



1 



(34)< 



2 \ e 



-S12 



2 \ e 



-S23 



2 / fi 



2 \ e 



e \-s23 



(12) (23) (45) 



Ls-if^,^) -2 



(3|15|4) 



L °Ni) 



(34) 



S45 



(IV.2) 



(IV.3) 



V sc = — 



1 / fi 



2 \ e 



2e \-s 23 



+ 1, 



(34)(3|15|4) ^0^—; 1 ( 3 |i5|4) 



\2 Li 



+ 



-S23 



(IV.4) 



(IV.5) 



(12) (23) (45) s 45 2 (12) (23) (45) sj 5 

Similarly, for the primitive amplitude A 6 (l+, 2^~ , 3 + , 4J , 5+), which contributes only at sub- 



leading order in N c , the tree amplitude is 



A 



tree 



(24)' 



(23) (31) (45) 
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(IV.6) 



The cut-constructible and scalar contributions are 



e 2 \ _s i2 J £ 



pec _ _ 



(IV.8) 



_JM_ Ls (s,, i (24) ((12) (34) -(14) (23)) ( ^ ^\ 

(23) (31) (45) H -»«'-»«>/ (23)(13) 2 (45) ^-'u' 

[13] (14) (24) L °(5ff) 



+ 2. 



(13) (45) S45 



vsc = i /j^_y + i (IV 9) 



2e V-S45 7 2 



(14) 2 (23\ /_ S12 x 1 (4|1|3) 2 (23) L i(^fS) , (14) 2 (2|3|1) L o(^f) 



(13) 3 (45) LS 1 1 J 2 (13) (45) s 2 3 + (13) 2 (45) s 23 

(2|1|3)(4|3|5) L i(^7t) (2|13|4) (14) L o(^ft) 1 [35] ([13] [25] + [23] [15]) 
(13) s\ 2 (13) 2 (45) s 12 2 [12] [23] (13) [45] 

(IV.10) 

The contributions with a closed fermion or scalar loop are rather simple. By Furry's theo- 
rem, the cases with vector-like couplings vanish, so A§ s, ^(l q , 2q, 3) = 0. The axial contribution 
y4| x (l+, 2^~, 3 + ) are also simple [49] and are given by 



F » = _ [ 53] [ 31](24)(^ip- i ^). (IV.ll) 
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